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^ ■ ABSTRACT 
P-f 

Q , We present the first search for galaxy counterparts of intervening high-z (2 < z < 3.6) sub-DLAs and DLAs towards GRBs. Our 

^ ■ final sample comprises of five intervening sub-DLAs and DLAs in four GRB fields. To identify candidate galaxy counterparts of 

the absorbers we use deep optical and near-infrared imaging, and low-, mid- and high-resolution spectroscopy acquired with 6 to 
[ 10-m class telescopes, the Hubble and the Spitzer space telescopes. Furthermore, we use the spectroscopic information and spectral- 

energy-distribution fitting techniques to study them in detail. Our main result is the detection and spectroscopic confirmation of the 
galaxy counterpart of the intervening DLA at z = 3.096 in the field of GRB 070721B (zgrb = 3.6298) as proposed by other authors. 
We also identify good candidates for the galaxy counterparts of the two strong Mgii absorbers at z = 0.6915 and 1.4288 towards 
^ I GRB 050820A (zgrb = 2.615). The properties of the detected DLA galaxy are typical for Lyman-break galaxies (LBGs) at similar 

OO . redshifts; a young, highly starforming galaxy that shows evidence for a galactic outflow. This supports the hypothesis that a DLA can 

OO ■ be the gaseous halo of an LBG. In addition, we report a redshift coincidence of diff'erent objects associated with metal lines in the 

' same field, separated by 130-161 kpc. The high detection rate of three correlated structures on a length scale as small as ~ 150 kpc 

I in two pairs of lines of sight is intriguing. The absorbers in each of these are most likely not part of the same gravitationally bound 

, structure. They more likely represent groups of galaxies. 

, Key words, galaxies: evolution - galaxies: formation - Galaxies: individual: DLA J0212-021 1 - quasars: absorption lines - quasars: 

CN ■ individual: QSO J1408-0346 - Gamma-ray burst: individual: GRBs 050730, 050820A, 050908, 070721B 

; 1. Introduction 

^ . Intervening absorption-line systems found in quasi- stellar-object 

■ Send offprint requests to: S. Schulze, sts30@hi.is (quasar, QSO) spectra play an important role in the observa- 
* Based in part on observations collected at the European tional study of galaxy formation and evolution. Unlike emission- 
Organisation for Astronomical Research in the Southern Hemisphere, selected objects, QSO absorption-line systems probe structures 
Chile, as part of the programs 075.A-0603, 075.A-0385, 077.A-0312, over the full mass range fro m under-dense regions to massiv e 
084.A-0303, 177.A-0591 and 275.D-5022. Based in part on obser- virialised dark matter haloes (lRauchlll998l: IWolfe et al.ll2005h . 
vations made with the NASA/ESA Hubble Space Telescope, ob- The most neutral-hydrogen-rich absorption systems are called 
tained at the Space Telescope Science Institute, which is operated damped Lyo^ absorbers (DLAs), if the neutral hydrogen col- 
by the Association of Universities for Research in Astronomy, Inc., density, A^(Hl), exceeds 2 X lO^^ cm'^ (log A^(Hl) > 
under NASA contract NAS 5-26555. Based in part on observations on o\ ^ u tat a / ^- t ^• u 

A -.u o V o rj. J u- u • . J u T . 20.3) and sub-DLAs (or sometimes super-Lyman-limit sys- 

made with the Spitzer Space Telescope, which is operated by the Jet x . ^ , , , . . . f i - ^ 

Propulsion Laboratory, California Institute of Technology, under a con- t^^^)/ ^^e column densi ty in neutral hyd rogen is between 

tract with NASA. Based in part on observations obtained at the Gemini ~ and 2 X 10^^ cm"^ (|Wolfe et al. ||2005', and ref. therem). 

Observatory, which is operated by the Association of Universities for DLAs represent the mam re servoir of neutral hydrogen m the 

Research in Astronomy, Inc., under a co-operative agreement with the high-z (z > 2) Universe (IWolfd Il986l: iPeroux et aH 120051: 

NSF on behalf of the Gemini partnership. IWolfe et alj l2Q05l) . Since the advent of QSO absorption-line 
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spectroscopy, o ver 1000 intervening sub-DLAs and DLAs have 
been detected (Prochaska et alJ l2005l: INoterdaeme et alJ 120091: 
Hipchaska & Wolfe 2009). 

Although their frequency and chemical composition are well 
known, the nature of their galaxy counterparts (hereafter called 
DLA galaxies) has remained poorly constrained for many years. 
Various models for DLA galaxies ex ist: e.g. rapidly-rotating 
proto-galactic disks at hig h reds hift CProchaska & Wolf3ll997l 
[l998: Wolfe & Prochaska 1 19981) . low surface brightness galax- 
ies ( Jimenez et al. 1999), faint and small gas-rich dwarf galaxies 
d^son 1988; Haehnelt et al. 1998; Okos hi & Nagashim a 2005), 
compact galaxies (Nagamine et al. 2007 i), and gaseous haloes of 
high-zLBGs (Fynbo et al. 1999; M0ller et al ."20021; |Fi^nbo_et^ 
l2008l:lRafels ki et al. 201 1). Chen & Lanzetta (2003 showed that 
the majority of low-z sub-DLAs and DLAs are late-type galax- 
ies and only a few are elliptical or dwarf galaxies. The conclu- 
sion that these low-z absorbers are indeed sub-DLAs or DLAs 
is drawn indirectly. Below z < 1.5, the hya absorption line is 
below 3000 A, a spectral rang^that is inaccessible with ground- 
based telescopes. lEllisonl (l2006l) showed that Mgii absorbers 
with £\yrest(Mg 11/12796) > 1 A, so-called strong Mgii ab- 
sorbers, are likely sub-DL As or DLAs (see also lRao & TurnsheS 
120001: lElhson et al.ll2009h . 

The successfully identified DLA galaxies allow the study of 
objects at the very faint end of the galaxy luminosity function 
(LF), objects that are usually missed in galaxy surveys. During 



carried out, but only with limited success (e. 


2. Smith et al. 


1989 


'Diorgovski et al."'l996>; "Le Brun et al."l997 


; 'Fvnbo et al.l 


1999 



also IBouche et al.l 120 12). Up to now, only about 17 sub-DLA 
and DLA galaxies have been identified below the redshift of 
unity (e.g. IPeroux et al.ll201 ll and reference therein) and more 
than 80 DLA ga laxy candidates have been reported to date (e.g. 
iRao et al. l l201lh . For the majority of these DLA galaxy candi- 
dates there is no spectroscopic confirmation yet. At high redshift 
the situation is worse; only 9 intervening sub-DLAs and DLAs 
with redshifts betwe en 2 and 3.15 have a confirmed ^^l^xy 
counterpart (high-z; Diorgovski et al. 1996; M0ller et al. 2002; 
Christensen et al. 2004b; M0ller et al. 2004; Weatherlev et al. 



2005; Fvnb o et al.|[2QM I^OTU Peroux et al. 2012; Bouche^etafl 
2012: .Krogager etal.ll2012t INoterdaeme et al...2012. and refer- 



ences therein). In addition to those, a few associate d DLAs, i.e. 
Zabs ^ Zem. h avc detected gala xy counterpar ts (e.g. Fvnb o et al 



1999 : lM0ller & Warren 1993i: lElHson et al.h2002: Miller et al 
2004l:lAdelberger et al.ll20Q6L 



Apart from the faintness of the galaxy counterparts, a fur- 
ther issue hampering the search for DLA galaxies is the glare 
of the bright background QSOs. Counterparts with a small an- 
gular distance from the QSO (impact parameter) are difficult 
to recover. This may lead to a possible misidentification of the 
DLA galaxy and hence an overestimation of its size, as DLAs 
at small impact parameters are easily missed. A diff'erent class 
of background light sources, gamma-ray bursts (GRBs), pro- 
vide a complementary approach to the problem. A GRB is a 
transient phenomenon, outshining the known Universe in y- 
rays for a fraction of a s econd up to hundreds o f seconds (e.g. 
iKouveliotou et al.lll993l: IZhang & MeszarosI 12004) . This short- 
lived emission is followed by an afterglow that can usually be 
detected from the X-rays over optical/near-infrared (NIR) wave- 
lengths to the radio for several weeks (e.g. Racusin et al. 2009; 
iKann et al.|[2QlQl: IChandra & Frailll20T2 ). An afterglow can out- 
shine the brightest quasars by several orders of magnitudes, but 



only for a very short period of tini e (~ 1 hr 1 day; e.g. 

iBloom et al.l[2QQ9l: iKann et al.ll2010l) . Similar to QSOs, GRBs 
can be found over most of the observable Universe, the most 
distant spectrosco pically -confir med burst being GRB 090423, 
at z ~ 8.3 CTanvir et alJ l2009l: ISalvaterra et al.ll2009l) . In fact, 
since most GRBs are associated with th e death of massive stars 
(Hiorth &Bloo^l2011l: IWooslevl l201lV this should allow us to 
detect them at higher redshifts than QSOs. In short, GRB after- 
glows can be brighter than QSOs but they are ephemeral, i.e. 
they vanish within a couple of months. This leaves the line- 
of- sig ht clear and without any interference from a bright object 
(e.g. Masetti et al."2003'; 'Vreeswiik et al.'"2003; Jakobsson et al.J 
2004: ElHson et al 2006: Prochter et al...2006: .Henriksen.2008i: 
Pollack et al.ll2009l) . 



GRBs are important for DLA studies for another rea- 
son. Nearly all GRB host galaxies have a DLA (GRB-DLA; 
iFvnbo etal.1 l2009l) . The distribution of their Hi column den- 
sity, and their metal-line strength distributions of Si ll /11526 and 
C IV /1/11548&1550 overlap with those of intervening DLAs, but 
also extend to larger values (iProchaska et al.ll2007bl:lFynbo et all 
2Q09). On average GRB -DLAs have larger metallicities than 
intervening DLAs (iFvnbo et al.l l2006l: iProchaska et al.l l2007bl: 
'Fynbo et al.l l2008l) . The diff'erence between in-situ and inter- 
vening DLAs lies in the way they prob e their host galax- 
ies t Vreeswiik et al. 20M IProchaska eTaLll2007bl: iFvnbo et al.l 
I2008L 2009). As GRBs are thought to originate from the col- 
lapse of a massive star, GRB -DLAs probe the line of sight to 
the location of a massive star, irrespective of the DLA orien- 
tation relative to the host galaxy if the progenitor is located in 
the DLA. Consequently, GRB -DLAs are selected by their star- 
formation rate (S FR). In contrast, the properties of an interven- 
ing DLA depend on the geometry of the DLA galaxy and its ori- 
entation relative to the line of sight. These are selected by their 
covering fraction and hence their cross-section, (t(Hi). A de- 
bated question in observational galaxy formation and evolution 
is the faint-end slope of the z ~ 3 galaxy LF. DLA galaxies probe 
the faint-end slope, but the in-situ and intervening DL A LFs are 
only identical if cr(lll) oc SFR dChen et al.l 12000: Fv nbo et al.l 
^008). The studies of both DLA populations therefore comple- 
ment each other well. 

Currently, seven intervening sub-DLAs and DLAs have 
been found in six GRB afterglow spectra, with absorber 
redshifts ranging between 2.077 and 3.564 (GR Bs 05073 0, 



050820A, 050908, 050 9 22C, 060607A, 07 0721B; Chen et al 
20081: 



2005 



Fox et al 



20091: IFvnbo et al.l l200^ 



Piranomonte et al.. ^ 
IVergani et al.l l200^ 



l20( 
I). 



^008; Chenetal 
Although the 



number of intervening (sub-)DLAs towards GRBs does not in- 
crease the number of known (sub-)DLAs significantly, the tran- 
sient nature of the background afterglow simplifies the search 
for their galaxy counterparts. After the afterglow has faded, deep 
follow-up observations are usually carried out to find the GRB 
host galaxy, typically reac hing a li miting magnit ude of ~ 2 7 
mag in the 7?-band (e.g. Hior th et al.i r2012: Malesani etani2012h . 
These observations are also suitable for the search for DLA 
galaxies. 

In this paper we present the findings from our search for the 
photometric counterparts of intervening sub-DLAs and DLAs in 
six GRB lines of sight. In Sect. [21 we introduce our methodology, 
present the sample selection, and describe how the data were 
reduced. The results are then presented in Sect. [3l confronted 
with current models of sub-DLAs and DLAs, and compared to 
low-z and high-z DLA galaxies in Sect.lH In Sect. [5] we draw our 
conclusions. 
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Table 1. Summary of photometric and spectroscopic data 
Photometry^ 



Instrument 


Filter 


Exposure 
time (s) 


Seeing 


Photo- 
metric? 


GRB 050730: RA = 


M^OS^^nni DEC : 


= -03°46' 17^^70 (J2000) 




VLT/FORS2 


R 8x 


250 + 24 X 300 


0^/6 


yes 


HST/ACS 


F115W 


6 X 1307 






VLT/ISAAC 


Ks 


32x60 


an 


yes 


Spitzer/IRAC 


3.6&5.8//m 


7200 






GRB 050820A: RA 


= 22^29^^38^14, DEC = 19°33'36'/61 (J2000) 




Keck/LRIS 


g 


2620 


Off! 


yes 


HST/ACS 


F625W 


2238 






VLT/FORS2 


R 


4x500 


0^/8 


no 


HST/ACS 


F115W 


4404 






HST/ACS 


F850LP 


14280 






Magellan/PANIC 


H 


8460 


0'/6 




VLT/ISAAC 


Ks 


32x60 


an 


yes 


Spitzer/IRAC 


3.6&5.8//m 


3600 






GRB 050908: RA = 


0lh21"^50^73, DEC 


= -12°57a7'/30 


(J2000) 




VLT/FORS2 


R 


12 X 500 


an 


yes 


VLT/ISAAC 


Ks 


32x60 


af6 


yes 


Spitzer/IRAC 


3.6&5.8//m 


7200 






GRB 070721B: RA 


= 02^12^^32^95, DEC = -02°ir4a'80 (J2000) 




VLT/FORS2 
HST/ACS 
Magellan/PANIC 
VLT/ISAAC 


R 

F115W 
H 

Ks 


24 X 235 
6 X 1307 
10320 
32x60 


af5 
a/5 

I'/O 


no 

yes 



Notes. Co-ordinates were taken from lMalesani et all (120121) . Observing 
conditions of Keck and Magellan data are quoted from IChen et al] 
12009). Column "Photometric?" states if observing conditions were 
photometric. Date of observations. GRB 050730: VLT/F0RS2 

- 2006/03/04 - 05/25, HST/ACS - 2010/06/10, VLT/ISAAC - 
2006/03/20, Spitzer/IRAC: 2008/03/11; GRB 050820A: Keck/LRIS 

- July 2006, HST/ACS - F625W, F115W, F850LP: 2005/09/26 
and 2006/06/05-11, Magellan/PANIC: August 2007, VLT/F0RS2 - 
2006/05/24, VLT/ISAAC - 2006/05/20, Spitzer: 2007/12/23; GRB 
050908: VLT/F0RS2 - 2007/07/21-08/14, VLT/ISAAC - 2007/07/11, 
Spitzer/IRAC: 2008/08/18; GRB 070721B: VLT/F0RS2 - 2007/10/04- 
11/13, HST/ACS - 2010/11/13, Magellan/PANIC - August 2007, 
VLT/ISAAC - 2007/09/22. 



Throughout t he pap er we refer to the Solar abundance com- 
piled in Aspl und et alJ (12009) and adopt cm"^ as the linear unit 
of column densities, N. We assume a ACDM cosmology with 
Ho = 71kms"^Mpc"^Qm = 0.27, and Qa = 0.73 (.Larson et al.i 

Hani). 



2. Data gathering and reduction 

2.1. Sample selection and data gathering 

All six GRBs with intervening sub-DLAs and DLAs were a part 
of deep imaging campaigns dedicated to detect their host galax- 
ies. Among these, we selected those that had sufficient spectro- 
scopic data, e.g. spectra obtained with different position angles 
(PAs), or multi-filter data to construct the spectral energy distri- 
bution (SED) of candidates and determine their nature. This re- 
duced the set to GRBs 050730, 050820A, 050908 and 070721B. 

In the following, we briefly describe the photometric and 
spectroscopic data analysed, as summarised in Table [T] The ob- 
serving cond itions durin g Keck and Magellan observations are 
described in IChen et al.l (2009), while those of other ground- 
based observations are summarised in Table [T] The astrometric 



uncertainty is ~ (y.'3; the uncertainty between the optical and 
NIR astrometry is similar, allowing us to unambiguously iden- 
tify objects in the diff'erent bands 

2.1.1. GRB 050730 

Afterglow spectra of GRB 050730 were acquired with 
VLT/F0RS2 with the 300V grating and with VLT/UVES 
with the red and blue arm centred at 3460 A and 5800 A, 
respectively (see Table [T]). In addition, a further spectrum was 
acquired with VLT/FORSl (grating 600V; RA=22.5?), as a part 
of the TOUGH (T he Optically Unbiased GRB Host) survey 
(lHiorthetal.112012) and three further spectra with VLT/F0RS2 
(PI: Ellison; Table [T]). For the latter, the slits were centred on the 
afterglow position, and spectra with difl'erent slit orientations 
were obtained with the 600B grism. In doing so, most of the 
field of view close to the GRB position is covered including 
several of the brightest close-by objects. This strategy enables 
us to localise emiss i on line object s via triangulati on, similar t o 
"M0ller et al. (2004), Fvnbo et al.' (^010), Fvnbo etaD (l201lh . 
.Krogager et al. (2012) Noterdaeme et al. (2012)B 

iFynbo et al.ll2QQ9l) reported the serendipitous discovery of 
a quasar at the same redshift as the intervening sub-DLA to- 
wards GRB 050730, 1775 south of the afterglow position. The 
quasar is blended with a 18.2-mag bright foreground star. In 
order to investigate the QSO- sub-DLA connection and to con- 
strain the quasar radiation field, we acquired a spectrum of the 
quasar with VLT/X-shooter (Table [T]). The field was also the 
target of deep imaging campaigns with VLT/F0RS2 (7?-band) 
and VLT/ISAAC ( T^s-band) as a part of the TOUGH survey 
(lHiorthetal.1 12012) . HST/ACS (F775W-band, PI: Levan) and 
Spitzer/IRAC (3.6 and 5.8 fim; Table[I]). 

2.1.2. GRB050820A 

Similar to GRB 050730, GRB 050820A was a target of an ex- 
tensive spectroscopic campaign. After the afterglow faded, sev- 
eral spectra were acquired with VLT/F0RS2 by centring the slit 
on the afterglow position and using three diff'erent slit orienta- 
tions, each one observed with two diff'erent gratings (300V and 
1028z; PI: Ellison), see Table [T] Slits were oriented to cover 
most of the area around the afterglow position and several of 
the brightest objects close to the afterglow position. The total in- 
tegration time of the data acquired with the 1028z grism is only 
half of the time spent on 300V, and the 1028z spectral range 
does not extend the range of the 300V data substantially, but its 
resolution is much higher than that of the 300V data. Hence, 
we do not use the 1028z data in our analysis. In addition, the 
field was observed in nine filters extending from 473.1 nm (g' 
band) to 5.8yL/m with VLT/F0RS2 and VLT/ISAAC as a part of 
the TOUGH survey, HST/ACS, KECK/LRIS, Magellan/PANIC, 
and Spitzer/IRAC (Table[B. 

2.1.3. GRB 050908 

For this burst, the amount of available data is small. An afterglow 
spectrum was obtained with the 300V grating (slit width of VfO 
and PA = 0^0) utiHsing VLT/FORSl, and with Gemini/GMOS- 
N UtiHsing the B600 grating (slit width of a.'75 and PA = 1 10^0; 
Table [B. In addition to these data, we included a spectrum ob- 
tained with VLT/FORSl (600B) and deep R- and /^s-band im- 



^ In fact, the galaxy counterparts of most sub-DLAs and DLAs at 
z> 2 were detected via triangulation. 
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Spectroscopy^ 



Instrument 


Grating 


Spectral 
range (A) 


Resolving 
power 


Exposure 
time (s) 


Slit 
width 


Position 
angle 


Seeing 


Photo- 
metric? 


GRB 050730 


VLT/FORSl 


600V+GG435 


4430-7370 


990 


5 X 1316 


r/3 


22?5 


a' 8 


yes 




600B 


3300-6210 


780 


4 X 1316 


r/3 


-20?0 


0'/9 


yes 


VLT/FORS2 


300V 


3400-9500 


440 


600 + 1200 


r/0 


0?0 


r/5 


no 


600B 


3300-6210 


780 


5 X 1316 


r/3 


33?5 


r/2 


yes 




600B 


3300-6210 


780 


4 X 1316 


r/3 


95?0 


r/3 


yes 


VLT/UVES 


346 and 580 
437 and 860 


3030-6840^ 
3730-10600^ 


~ 40000 


3000 


r/o 


0?0 


r/8 
r/1 


no 
yes 




UVB 


3000-5500 


4000 


2x900 


r/3 








VLT/X-shooter^ 


VIS 


5500-10000 


6700 


2x900 


r/2 


43?7 


a' 8 


yes 




NIR 


10000-25000 


4000 


6x300 


r/2 








GRB 050820A 




300V 






4 X 1316 


r/3 


-52?0 


r/2 


yes 




300V 


3400-9500 


440 


8 X 1316 


r/3 


34?0 


r/3 


yes 


VLT/FORS2 


300V 
1028Z+OG590 






4 X 1316 
2 X 1316 


r/3 
r/3 


83?0 
-52?0 


0^/7 
0^/8 


yes 
yes 




1028Z+OG590 


7730-9480 


2560 


4 X 1316 


r/3 


34?0 


r/2 


yes 




1028Z+OG590 






2 X 1316 


r/3 


83?0 


0^/7 


yes 


VLT/UVES 


390 and 564 
437 and 860 


3260-6680^ 
3730-10600^ 


~ 40000 


1800+1815 
2430 


r/o 


0?0 


r/o 


yes 


GRB 050908 


Gemini/GMOS-N 


B600+G0305 


5118-7882 


844 


2 X 1200 


^'75 


110?0 






VLT/FORSl 


300V 


3400-9500 


440 


3 X 1200 


r/o 


0?0 


^'6 


yes 


600B 


3300-6210 


780 


6 X 1345 


r/3 


53? 8 


r/2 


yes 


GRB 070721B 


VLT/FORSl 


300V 


3400-9500 


440 


5 X 2735 


r/3 


-118?5 


r/2 


yes 


VLT/FORS2 


300V 


3400-9500 


440 


4x 1800 


r/o 


-118?3 


r/2 


yes 



Notes. O bserving conditions of the afterglow spe ctros copy campaigns, and the VLT/UVES data of GRB 050730 and of GRB 050820A are 
taken fromlF vnbo etaP (120091) . iLedoux et al.1 (|2g09l), andl Vergani etall (120091). The VLT/UVES and VLT/X-shooter spectrographs consist of two 
and three arms, respectively, displayed in the column "Grating". Observations in the different arms are carried out simultaneously. Resolving 
power of VLT/FORSl +2 and Gemini grisms are given for a slit width of V . The stated resolving power of the X- shooter NIR spectrum was 
decreased during the data reduction to increase S/N. The effective resolving power of the X- shooter NIR arm is around ten times lower. Date 
of observations. GRB 050730: VLT/FORSl 600V+GG435 - 2007/02/23-24, VLT/F0RS2 600B - PA = 2(f.O: 2006/06/22, PA = 0^0: 2005/07/31, 
PA = 33?5: 2006/04/28-05/28, PA = 95?0: 2006/05/28, VLT/UVES: 2005/07/31, VLT/X-shooter: 2010/04/18; GRB 050820A: VLT/F0RS2 300V 
- PA = -52?0: 2006/07/22-08/16, PA = 34?0: 2006/08/31-09/22, PA = 83?0: 2006/08/17-09/15, 1028z - PA = -52?0: 2006/08/16, PA = 34°0: 
2006/09/17-22, PA = 83?0: 2006/08/17, VLT/UVES: 2005/08/20; GRB 050908: Gemini/GMOS-N - 2005/09/08, VLT/FORSl - 300V: 2005/09/08, 
600B: 2007/08/16-09/06; GRB 070721B: VLT/FORSl - 2007/11/16-12/04, VLT/F0RS2 - 2007/07/22. Spectral range of the blue and red arms 
are 3030-3880 A and 4760-6840 A, respectively. Spectral range of the blue and red arms are 3730-4990 A and 6600-10600 A, respectively. 

The target of the X-shooter observation was QSO J1408-0346. Spectral range of the blue and red arms are 3260-4450 A and 4580-6680 A, 
respectively. 

Tab.[T] — continued 



ages acquired with VLT/F0RS2 and VLT/ISAAC as a part of 
the TOUGH survey, respectively, and deep Spitzer/IRAC images 
acquired at 3.6 and 5.8 yum (Table [T]). 



2.1.4. GRB 070721 B 



We use the afterglow spectrum (PA=-118?3) acquired with 
VLT/F0RS2 (Table [B, and a spectrum obtained with 
VLT/FORSl (RA=-118?5) several months after the GRB faded 
as a part of the TOUGH campaign. In addition, several deep 
imaging campaigns targeted this field. We use the VLT/F0RS2 
7?-band and VLT/ISAAC K^-hmid data acquired as a part 
of the TOUGH campaign, HST/ACS F775W-band data, and 
Magellan/PANIC if -band data. 



2.2. Data reduction 
2.2.1. Imaging 

Ground-based data - VLT/FORS and Keck/LRIS data were 
reduced in a standard way, VLT/FORS data with I RAF and 
Keck/LRIS data with IDL (Todv 1986: see lMalesani etal.ll2012 
for rn ore details on the VLT data reduction, and iPerley et al. 
[2012' on the Keck data reduction). The NIR data acquired with 
VLT/ISAAC were red uced in a standa rd way with the jitter 
routine in eclipse (Dev illardI Il997|). The red uction of the 
Magellan/PANIC data is presented in lChen et al.] (120091 

HST data - HST/ACS images of GRB 050730 and 070721B 
consist of 6 dithered exposures in the F115W filter. Individual 
exposures, after standard "on-the-fly" processing, were retrieved 
from the STScI archiveQ These were subsequently cleaned for 
bias striping, introduced due to the replacement electronics after 



http://archive.stsci.edu 
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Servicing Mission 4 (May 1 1-24, 2009), and then drizzled using 
the multi drizzle software into final science images. At this 
stage we adopt a pixel scale of 0^/033, approximately 2/3 of the 
native pixel scale. T he data redu ction of the GRB 050820A HST 
data is presented in et al.l (I2OO9.) . The pixel scale of these 
drizzled images is O'.'OS in contrast to that of the images of GRBs 
050730 and 070721B. 

Spltzer Data GRBs 050730, 050820A, and 050908 were ob- 
served with the I nfrared Array Camera aboard Spitzer (IRAC; 
iFazio et aD 12004 ) in the bandpasses centred at 3.6 and 
5.8 )um. We downloaded the processed PBCD data from the 
Spitzer Heritage Archive and followed the basic procedures for 
aperture photometry in the IRAC Instrument Handbook0 using 
an extraction radius of 2 native pixels (2V4, or 4 pixels in the 
default archive resampling). 

Photometry - In order to measure the total flux, we 
applied aperture phot ometry, utilising Source Extractor 
(iBertin & ArnoutsI Il996h . and applied an aperture correction 
assuming a stellar point-spread fu nction (PSF; for th e aperture 
correction of the HST images see S irianni et aD [2QQ5b . 

The quasar in the field of GRB 050730 is strongly blended 
with a foreground star at a projected distance of V!2. To perform 
reliable photometry, the contaminating star was subtracted from 
the VLT/FORSl images. A model of the FORSl PSF was con- 
structed from ~60 field stars using Daophot as implemented in 
I RAF. This PSF model was then fitted simultaneously to the star 
and the QSO. To test the quality of the fit, we subsequently sub- 
tracted the stellar contribution from the input image. Photometry 
was then performed in the final star- subtracted image. 

Instrumental magnitudes obtained from optical ground- 
based data_were_£alibrat^ 2-3 photometric standard 
stars (lLandoltlll99i IStetsonl bOOO). NIR magnitudes were cal- 
ibrated with more than three 2MASS stars and HST magnitudes 
against zeropoints, computed from their FITS headers. 

2.2.2. Spectroscopy 

VLT/F0RS1 and F0RS2, Gemini/GMOS and VLT/UVES data 
- FORS and Gemini data were reduced in a s tandard way with 
I RAF (for more details see iFvnbo et al.ll2009l). The data reduc- 
tion of the UVES data is described in lLedoux et al.l (l2009l) . 

VLT/X-shooter data - The QSO in the field of GRB 050730 
was observed with X-shooter. These data were reduced with the 
X-shooter pipeline vl.2.2Q The 1-dimensional spectrum of the 
QSO {R = 20.76 mag; corrected for Galactic extinction) in the 
field of GRB 050730 could not be extracted with the X-shooter 
pipeline, because of being blended with the 18.2-mag bright 
foreground star Vf2 SW mentioned above (see Fig.[T]). 

To estimate the flux of the quasar over the entire wavelength 
range of the X-shooter spectrograph, we proceeded in the fol- 
lowing way. For each arm, we first extracted a profile of the su- 
perposed spectral PSFs at several wavelengths. Each profile was 
best described by two Gaussians, representing the QSO and the 
star. The centre of their peaks and the widths did not change 
with wavelength, allowing us to fix these parameters to their 
mean values for each arm of the X-shooter spectrograph. Then, 



^ http :// http : //sha.ipac .caltech.edu/applicat ions/Spitzer/SHA/| 
^ http ://w w w. eso . org/sci/soft ware/pipelines/ 



we extracted the profile for every wavelength and fi t it with the 
aforementioned model, using the routine mpfit by Markward^ 
( 2009), in IDL. To identify outliers, a 9-point median profile was 
fitted for every wavelength, i.e. four points red ward and four 
points blueward of the considered wavelength. Data points were 
rejected if they deviated by more than 3cr from the median spec- 
tral PSF of the quasar and the star. 

The QSO was barely visible in the NIR spectrum. We there- 
fore determined the centre of peak and the FHWM of the 
star and set the QSO centre of peak to its expected position. 
Furthermore, we reduced the spectral resolution to increase the 
signal-to-noise ratio (S/N) by rebinning the spectrum. The un- 
certainties in the FWHMs and the centre of peaks varied be- 
tween 0.06 and 0.39 px and 0.3 and 0.18 px for each arm (pixel 
scale: UVB and VIS Qf:i6/px, NIR a'21/px), respectively, smafl 
enough to be neglected. 

We compared the quality of our method with the MCS 
decon yolution technique develop ed by Magain, Cou rbin, & 
Sohv (iMagain et al.iri998h used in lLetawe et al.l (l2008l see also 
ICourbin et al.ll200Ql) . This method requires an unblended star to 
be observed with the identical instrument setup. Because of the 
short slit length, no such object was covered by the slit. For a 
rough estimate, we used the observed standard star as reference, 
although it was observed with a slit width of 5''. Within the er- 
rors both methods give the same results. 

To flux calibrate the spectrum, we reduced the data of the 
standard star GD50 with the X-shooter pipeline. The ID spec- 
trum was extracted using the routine apall in optimal ex- 
traction mode. The spectrum was then divided by the corre- 
spon ding reference spectr um from the CALSPEC HST database 
(Bo hlin & Gillilandl [20041) to deduce the response function. In 
addition, we corrected the flux-calibrated standard star spectrum 
for undulations by smoothing the ratio between the observed 
flux-calibrated and expected spectra with a Hamming filter (win- 
dow size: 40 A in the UVB and VIS arm) in regions that were 
not aff'ected by strong telluric lines or stellar absorption features. 
We then applied the corrected response function to the QSO and 
used the acquisition i mage to secure the ab solute flux calibration. 
Finally, we followed lCardeUi et al.1 (1 19891) to correct the QSO for 
Galactic dust attenuation (E(B - V) = 0.048 mag). No attempt 
was made to correct for telluric absorption lines due to the lack 
of a suitable telluric standard star observed the same night. This 
has no implications on our analysis. 

2.3. SED fitting 

The photon ietric SEDs of the galaxies wer e modeled within 
Le Phare (lArnouts et al.1 1 19991: lllbert et al.l j2006)PI Here, we 
used a grid of galaxy templates based on feruzual & CharlotI 
( 2OO 3I) stellar p opulation synthesis mo dels with a Chabrier IMF 
(Chabrieret al. 2000) and a Calzetti (Calzetti et al.l l2000h dust 
attenuation curve. In those cases, where prior information about 
the redshift was available from spectroscopy, we fixed this pa- 
rameter to the respective value. For a description of the galaxy 
templates, physical paramet ers of the gal axy fitting and their 
error estimation we refer to iKrtihler etal.! (120111) . To account 
for zeropoint off'sets in the cross calibration and absolute flux 
scale, an systematic error contribution of 0.05 mag was added in 
quadrature to the uncertainty introduced by photon noise. 



^ http://www.cfht.hawaii.edu/ arnouts/LEPHARE 
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Table 2. Properties of the intervening sub-DLAs and DLAs 



GRB 


^GRB 


-^abs 


IVJg i V 1 1^ 


[Si/H] 


(A) 


Detected lines 


References 






3.56439 


20.3 + 0.1 


< -1.3 


< 1.01^ 


Al II, Fe II, Si II, Si ii* , C iv. Si iv 




050730 


3.969 


3.02209 


19.9 + 0.1 


-1.5 + 0.2^ 




Al II, Fe II, Si II 


1,2,3 


050820A 


2.615 


2.3598^ 


20.1+0.2 


-1.5 + 0.2 


0.17 + 0.01^ 


Fe II, Si II, C IV, Mg ii 


4 


050908 


3.3467 


2.6208 


20.8 + 0.1^'^ 


> -1.40^ 


2.24 + 0.06 


Al II, Fe II, Si II, C IV 


2 


070721B 


3.6298 


3.0939^ 


20.1 +0.3^'^ 


> -0.66^ 


2.35 + 0.69^ 


Al II, Fe II, Si II, C IV 


2 



Notes. This line is blended. The EW listed is the value of the blend. Without ionisation correction. IVergani et al.l (I2OO90 reported 
£ Wrest (Mg II i2796) < 0.42 A with a m ost likely value of 0.3 1 A. This \york. iFvnbo et al.l (12009) suggested that the absorber is a sub-DLA. 

Derived in the optically thin limit. IChen et al.l (120091) and lFvnbo et al.l (120091) proposed that a galaxy 0'.'9 SE of the afterglow position is the 
galaxy counterpart to the intervening absorber. 

References. fP lChen et al.l(l2005l) : (2) lFvnbo et an(l2009l) : (3) IStarling etan(l2005h : (4) IVergani et'an(l2009h 



Table 3. Properties of selected Mg ii absorbers 



GRB 


^GRB 


-^abs 


Wrest (MgiU2796) 
(A) 


Detected lines 


References 


050730 


3.969 


2.25313 
1.773F 


< 0.78 (0.65)"^ 
0.93 + 0.03^ 


Fe II, Mg I, Mg II 
Fe II, Mg I, Mg II 


1,2 
1,2 


050820A 


2.6147 


1.4288 
0.6915 


1.32 + 0.02 
2.87 + 0.01"^ 


Fe II, Mg I, Mg II 
Ca II, Fe II, Mg ii 


1 
1 


050908 


3.3467 


1.5481 


0.82 + 0.07"^ 


Fe II, Mg II 


3 



Notes. We only summarize the properties of those Mg ii absorbers that most likely have an intervening sub-DLA or DLA. Following the typical 
approach in the literature, we call absorbers with £ Wrest (Mgii/12796) > l.OA strong, with 0.3 A < EWrest (Mg ii /12796) < l.OA intermediate 
and with £Wrest (Mg II/12796) < 0.3 A weak systems. Absorbers are marked by ^ if the doublet is saturated. IChen et aLl (l2005h identified two 
velocity components that are separated by 57 km s~^ 

References. fP lVergani etaP (120091) : (2) lProchaska et al.l(l2007al) : (3) lFvnbo et all (120091) 



3. Results 

Afterglow spectroscopy allowed us to identify seven interven- 
ing sub-DLAs and DLAs in six GRB fields. For four of these 
fields, we had sufl&cient data to elucidate the nature of several 
galaxy counterpart candidates. In the following, we present how 
candidates were selected, and present the properties of the inter- 
vening sub-DLAs and DLAs towards GRBs 050730, 050820A, 
050908 and 07072 IB and our findings on the most likely galaxy 
counterpart for each field. A summary of the properties of the 
intervening sub-DLAs and DLAs is presented in Tables [2] and [3] 
Most of the shown values are taken from the literature. 

In order to select candidates that are close enough to the 
GRB line of sight to produce an intervening absorption-line- 
system in the afterglow spectrum, we had to set an upper limit 
on the extent of the possible DLA galaxies. Theoretically, it 
is very difl&cult to set a meaningful upper boundary, because 
diflTerent models exist for their galaxy counterparts. The ob- 
served impact parameter distribution of confirmed DLA galax- 
ies, based on Peroux et al. ( 201 1) and Krogager et al. ( 2012), ex- 
tends from 0.4 to 182 kpc and has a mean value of ~ 25 kpc. 
The impact param eter distribution of DLA galaxy candidates by 
iRao et al.l (l201ll) shows similar characteristics. The majority of 
DLA galaxies have a small impact parameter, however there are 
few cases with large impact parameters 100 kpc ). We there- 
fore follow the statistical approach by iRao et al.l (12011)) who 
used the galaxy number density as a function of impact param- 



eter as a criterion. They found that the galaxy number density is 
comparable to the number density of foreground and background 
galaxies, i.e. a chance association is more likely, if the impact 
parameter exceeds 100 kpc. As a first assumption, we limit our 
study to those candidates within 100 kpc from the GRB line of 
sight, keeping in mind that this value was derived for intervening 
DLAs between z = 0.5 and 0.8. 

3.1. GRB 050730 

GRB 050730 occurred at a redshift of z = 3.969. Its after- 
glow spectrum contains an intervening sub-DLA (log A/^ (H l) = 
19.9 ±0.1, [Si/H] = -1.5 ± 0.2)0 and an intervening DLA 
(logA^(Hi) = 20.3 ± 0.1, [Si/H] < -1.3) at z = 3.02209 and 
z = 3.56439, respectively (Tables [2l [3]). In addition the after- 
glow light traversed an intermediate and a strong Mg ii absorber 
at z = 2.25313 and 1.7731, respectively. 

In Figure [TJ we show the field of GRB 050730 and zoom-ins 
to the inner 28V 8 and 6V0; the region at which the impact pa- 
rameter exceeds 100 kpc is highlighted. At the redshift of the 
aforementioned absorption-line systems, the maximum impact 
parameter of 100 kpc translates to an transverse distance be- 
tween 12V1 and 13V9; the smaller value belongs to z = 1.7731 



^ The metallicity is not corrected for ionisation effects. 
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Table 4. Galaxy counterpart candidates of the intervening absorption line systems 



GRB 050730 (ORB: z = 3.969, DLA: z = 3.56439, sub-DLA: z = 3.02209, Mgii absorbers: z = 2.25313, 1.7731) 

Candidate Remark 6i") (mag) F775W (mag) Ks (msig) 3.6//m(mag) 5.8/im(mag) 



Al star 12.8 23.46 ± 0.03 23.08 ± 0.01 22.51 +0.25 23.30 + 0.20 

A2 12.3 28.50 ± 0.26 

A3 10.6 26.53 ±0.12 26.60 ± 0.09 

A4 star 7.3 23.88 ± 0.03 23.71 ± 0.02 

A5a 2.5 27.20 ±0.27 27.98 ±0.18 

A5b^ 2.1 28.16 ±0.19 

A5c^ 2.0 28.43 ± 0.22 

A5d 1.8 28.84 ± 0.30 

A6a 1.0 28.64 ± 0.26 

A6b 0.7 28.37 ± 0.24 

A6c 0.4 28.11 ±0.19 

A6d 1.4 28.64 ± 0.26 

A6e 0.9 28.64 ± 0.26 

A6f 1.5 28.29 ±0.20 

A7 star 4.2 23.38 ± 0.03 22.84 ± 0.01 22.28 ±0.15 22.78 ±0.15 

A8 6.3 27.49 ±0.12 

A9 9.8 27.38 ±0.12 

AlO 11.5 27.10 ±0.13 

All^ z~0.16?;z<2.81 13.3 

QSO^ z = 3.022 17.5 20.87 ± 0.05 -20.58 

Star^ star 18.3 18.20 ±0.04 



GRB 050820A (GRB: z = 2.615, DLA: z = 2.3598, Mgii absorbers: z = 1.4288, 0.6915) 



Candidate 


Remark 


en 


g (mag) 


F625W (mag) 


R (mag) 


F115W (mag) 


F850LP (mag) 


H (mag) 


Ks (mag) 


Bl^ 




11.8 




25.18 ±0.04 




24.87 ± 0.03 


24.90 ± 0.06 






B2 


z = 0.693 


12.3 


25.41 ± 0.07 


23.59 ± 0.03 


23.20 ± 0.03 


22.35 ± 0.01 


21.81 ±0.05 


20.81 ±0.10 


20.47 ± 0.05 


B3 


Star 


5.2 


25.24 ± 0.06 


23.90 ± 0.04 


23.65 ± 0.04 


23.01 ± 0.02 


22.57 ± 0.05 


21.41 ±0.14 


21.05 ±0.07 


B4 


2.8 


23.39 ± 0.05 


22.11 ±0.01 


21.97 ±0.03 


21.68 ±0.01 


21.34 ±0.01 


20.94 ±0.11 


21.19 ±0.07 


B5 


z = 0.9 ±0.1 


3.7 


26.68 ±0.14 


26.36 ± 0.09 


26.48 ± 0.26 


25.58 ± 0.06 


25.36 ± 0.06 






B6 


7 - 1 5+0-2 

z = 2.615 


4.9 


25.23 ± 0.06 


24.94 ± 0.09 


24.80 ± 0.09 


24.32 ± 0.04 


23.67 ± 0.05 


22.28 ± 0.21 


22.20 ±0.15 


B7« 


1.6 


25.75 ± 0.08 


25.99 ± 0.07 


25.66 ±0.10 


25.90 ± 0.06 


25.83 ± 0.08 






B8N« 
B8S« 


z = 2.615 


0.3 
0.1 




26.16 ±0.08 
26.19 ±0.08 




25.83 ± 0.08 
25.72 ± 0.07 


26.18 ±0.06 
25.82 ± 0.05 






B9 


z = 0.428 


3.4 


25.49 ± 0.07 


24.94 ± 0.05 


24.45 ± 0.07 


24.76 ± 0.03 


24.57 ± 0.04 






BIO 


star 


3.6 


22.11 ±0.05 


20.58 ± 0.00 


20.31 ±0.03 


19.69 + 0.00 


19.16 ±0.00 


18.60 ±0.05 


18.81 ±0.04 



Candidate 


3.6 //m (mag) 


5.8 yum (mag) 


Bl^ 


23.45 ± 0.30 




B2 


20.59 ± 0.05 


21.40 ±0.20 


B3 


20.94 ± 0.06 


21.04 ±0.18 


B4 


21.65 ±0.10 


21.68 ±0.32 


B5 






B6 


21.81 ±0.10 




B7« 






B8N^ 






B8S^ 






B9 






BIO 


19.51 ±0.03 


20.29 ± 0.09 



Notes. Observed magnitudes are given in the AB system (not corrected for Galactic extinction). VLT and Magellan Vega-magnitudes were 
converted into the AB system by adding: VLT//? 0.21 mag (Bla nton & Row eis 2007), YLT/K, 1.895 mag (ESQ), and Magellan/// 1.34 mag. The 
selective Galactic extinctions, E(B - V), are 0.051 and 0.044 mag, for GRBs 050730 and 050820A, respectively. Impact parameters were derived 
from the HST images, if available, otherwise from R-bmid images. The uncertainty in the impact parameter is O'.'S - 0'.'4. Spectra were extracted for 
boldfaced objects. Based on the spectrum, the SED or the shape and colour of an object, we report in the column "remark" if it is a star or at which 
redshift the galaxy is. The 3cr limiting magnitudes are: GRB 050730: R = 27.1, F715W = 28.6, K, = 22.8, 3.5//m = 23.5, 5.8yum = 21.5 mag, 
GRB 050820A: = 27.6, F625W = 27.8, R = 26.1, F115W = 28.1, F850LP = 28.0, H = 23.0, = 23.0, 3.5yL/m = 23.5, 5.8yL/m = 21.7mag, 
Blended object. The magnitude is only an estimate. The spectrum can be the integrated spectrum of the blended objects. The object consists of 
at least two galaxies and is detected in all filters. We do not report any magnitude of the individual objects or the compound in this work. We detect 
the continuum down to 4630 A. If this break is caused by Lya absorption, the redshift is z = 2.81, or if this is the redshifted 4000 A feature, the 
redshift is z = 0.16. The /?- band brightness of the quasar and the star were obtained from an afterglow image taken on 23 Fe bruary 2007 using 
PSF photometry (Sect. l2.2TI ). The blended star is saturated in the HST image. The PSF was subtracted with Tiny Tim ('Krist'"1993) assuming a 
FWHM of 5". While the wings of the star were partly removed, the residuals in the core are quite strong. The measurement of the QSO magnitude 
is only an estimate. ^"^^ The compound is a blend of several individual objects. We only report the brightness in the HST images of the brightest 
object. 



and the larger to z = 3.56439, including the astrometric uncer- The GRB was a target of an extensive spectroscopic cam- 
tainty of the afterglow localisation of (7/3. paign; several low-resolution spectra with a total of five diff'er- 

ent position angles (PA) were obtained with VLT/FORSl and 
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GRB 050908 (DLA: z = 3.3467, DLA: z = 2.6208, Mgii absorber: z = 1.5481) 



Candidate 


Remark 


en 


R (mag) 


Ks (mag) 


3.6 //m (mag) 


5.8 yum (mag) 


CI 


z ~ 2.71 


12.6 


27.10 ± 0.20 








C2 




11.5 


24.96 ± 0.04 




23.25 ±0.19 




C3 




6.6 


25.97 ± 0.09 








C4^ 




12.4 


26.17 ±0.09 








C5^'^ 




3.7 


25.05 ± 0.04 




23.0J + 0.23 


OO A A 1 f\ A A 

2ZA4 + (J. 44 


C6^'^ 




2.8 


25.39 ± 0.05 








C7^ 




1.1 


26.98 + 0.18 








C8 




5.5 


26.59 ± 0.21 








C9 




2.5 


25.96 + 0.12 




~ ZJ.U ± U.J 




GRB070721B (GRB : z 


= 3.6298, DLA: z = 3.0939) 






Candidate 


Remark 


en 


R (mag) 


F115W (mag) 


H (mag) 


Ks (mag) 


Dl 




12.8 




27.40 ±0.12 








z ^ 2.64 


3.0 




27.34 ±0.12 






D2b^ 


2.8 




28.23 ± 0.22 






D3 


z = 3.096 


1.0 


24.48 ± 0.02 


24.41 ± 0.02 


23.40 ± 0.07 


23.56 ± 0.23 


D4 


z = 3.631 


0.1 




27.31 ±0.15 






D5 




6.9 


26.82 ±0.18 


26.44 ±0.10 






D6a^ 


z ^ 2.64 


11.7 




25.05 ± 0.04 






D6b^ 


12.2 




26.90 ± 0.09 






D7 


z = 3.615 


20.7 


23.77 ± 0.01 


23.74 ± 0.03 


24.02 ±0.11 


23.20 ±0.17 



Notes. The selective Galactic extinctions, E{B - V), are 0.025 and 0.031 mag for GRBs 050908 and 070721B, respectively. Spectra were extracted 
for boldfaced objects. Based on the spectrum or the SED of an object, we show in the column "remark" if it is a star or at which redshift the 
galaxy is. The 3cr limiting magnitudes are: GRB 050908: R = 27.4, = 23.0, 3.6 /urn = 23.8, 5.8 yum = 22.3 mag and GRB 070721B: R = 27.5, 
F775W = 28.7, H = 23.1, Ks = 23.2 mag. Blended object. The magnitude is only an estimate. The spectrum is the integrated spectrum of the 
blended objects. We do not attempt photometry on objects which are strongly blended in the Spitzer data in this work, with the exception of 
objects C5, C6, and C7, which are sufficiently closely blended that we report the (summed) photometry of all three objects as a group. 

Tab.m — continued 




Fig. 1. Left panel: The field of view of GRB 050730 (z = 3.969) observed in F115W with HST/ ACS. Middle panel: Zoom-in to 
the inner 28V 8. The diff'erent sHt orientations from Table [T] are over-plotted. Right panel: Zoom-in to the inner 6''; the panel was 
smoothed by a boxcar filter (width: 0'.'132 = 4 pixel). The afterglow position is marked by an 'x' in all panels. The 100-kpc radius 
(shown in the middle panel) represents the assumed maximum impact parameter of the absorber at z = 3.56. Objects are labelled if 
they lie within the circle and if they are covered by any of the slits. Their magnitudes are summarised in Tabled 



VLT/F0RS2 (Fig.m Table [B- Within the 100-kpc radius, 19 
objects were covered by a slit. In Table |4] we summarise their 
magnitudes in diff'erent filters and their angular distances from 
the afterglow position. Four of these objects were bright enough 
to allow the extraction of a spectrum. The VLT/FORS 1 spectra 
of objects Al and A4 are of very low S/N. We detect their con- 
tinua down to 4500 A, but we detect no absorption or emission 
lines. If these objects were the galaxy counterparts to the sub- 
DLA and DLA, we expect to see the onset of the Lya forest at 
~ 4890 A and ~ 5549 A, respectively, which we do not observe. 
The Lyof non-detection and the extension of the continuum to 
even shorter wavelengths, rules them out as the galaxy counter- 



parts. We detect the continuum of object A7 at very low S/N in 
the VLT/F0RS2 data, but no absorption or emission lines. The 
objects Al, A4 and A7 are likely late-type stars based on their 
colours and the fact that the morphology and size of their PSFs 
(FWHM = (y.'ll) do not diff'er from point sources. Object All 
lies at the edge of the 100-kpc radius. We detect its continuum 
down to 4630 A in the VLT/F0RS2 data and a drop in flux blue- 
ward of it. Assuming that this is the B aimer break, the redshift 
of the galaxy is ~ 0.16, not in agreement with any of the inter- 
vening absorbers. 

In conclusion, we do not identify any galaxy counterpart can- 
didate of the intervening sub-DLA and DLA down to a limit- 
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Fig. 2. Similar to Fig. [H but for GRB050820A (z = 2.615). 
In contrast to ORB 050730, the field of GRB 050820A was ob- 
served in 9 filters ranging from 473 nm to 5800 nm, in addition to 
extensive spectroscopy campaigns. This allowed us to construct 
the SED of several objects in the field of view. The width of the 
displayed field is 29'', twice the maximum impact parameter of 
100 kpc of the strong Mgii absorber at z = 0.6915, while the 
circle displays the maximum impact parameter of 100 kpc at the 
redshift of the intervening sub-DLA atz = 2.3598. Every object 
within 100 kpc is labelled either if it is either detected in at least 
five filters or covered by a slit. Objects B7 and BIO appear to be 
not in the slit, due to the smaller PSF of the HST in comparison 
to the VLT. 



ing magnitude of F775W = 25.7 (24.6) mag and of any of the 
intermediate Mgii absorbers down to a limiting magnitude of 
F775W = 26.5 (25.1) mag, assuming a maximum impact param- 
eter 50 (100) kpc. These limits were calculated by considering 
all objects within 50 (100) kpc but excluding those for which we 
elucidated the nature or those with a stellar PSF. 

iFynbo et aP (l2009l) reported the serendipitous discovery of a 
QSO at the redshift of z = 3.023, very similar to the redshift of 
the sub-DLA in the afterglow spectrum. The QSO is 17V 5 south 
of the afterglow position (see Fig.[T] Tabled, corresponding to 
a projected distance of 136.7 kpc at z = 3.02209, the redshift of 
the sub-DLA towards the GRB. It has a brightness of 20.76+0.05 
mag in the /?-band (corrected for Galactic extinction. Table |4]), 
but is blended with a 18.2-mag bright K-type star. Strictly speak- 
ing, the impact parameter exceeds the assumed maximum im- 
pact parameter of 100 kpc. Given the redshift coincidence with 
the sub-DLA in the afterglow spectrum we discuss this corre- 
lated structure in Sect. 14.41 

3.2. GRB050820A 

iLedoux et al] (l2005h and lVergani et al.l (|2009l) reported an inter- 
vening sub-DLA at z = 2.3598 with log (H i) = 20. 1 ± 0.2 and 
a metallicity of [Si/H] = -1.5 + 0.2 and two strong Mgii ab- 



1 1 1 1 1 1 1 |: 1 1 1 1 1 


1 n 1 




ya 








1 1 1 1 1 1 1 K 1 1 1 1 1 1 1 



4300 4400 4500 

Observed wavelength (A) 



Fig. 3. The GRB 050908 afterglow spectrum, obtained with 
VLT/FORSl, centred on the intervening hya absorption line 
(^j3LA = 2.6208; logA^(Hl) = 20.8 + 0.1). A neutral hydrogen 
column density fit to the damped hya line is shown with a solid 
line, while the shaded region indicates the Icr errors. 



Fig. 4. Similar to Fig.[T] but for GRB 050908 (z = 3.3467). The 
diff'erent slit orientations from Table [T] are over-plotted. The cir- 
cle radius of 12V 7 represents the assumed maximum impact pa- 
rameter of 100 kpc of the intervening sub-DLA at z = 2.6208. 
Objects are labelled if they lie within the circle and if they are 
covered by any of the slits. Their magnitudes are summarised in 
Tabled 



sorbers between z = 0.6915 and 1.6204 towards GRB 050820A 
(ZGRB = 2.615, TablesElIll) 

In Figure O we show the inner 29'' around the afterglow po- 
sition and highlight the region at which the impact parameter is 
100 kpc. At the redshift of the intervening absorbers, the maxi- 
mum impact parameters correspond to 12! 4 for z = 2.3598 and 
14V4 for z = 0.6915, including the uncertainty of 0'.'4 in the af- 
terglow localisation. 

After the afterglow faded, several spectra were acquired with 
VLT/F0RS2. In addition, the field was also the target of an ex- 
tensive photometric campaign covering nine filters from 473.1 
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nm (^'-band) to 5.8 fim; for full details we refer to Sect. 12. H and 
Table [TJ This allows us to elucidate the nature of several objects 
within 100 kpc, not only of those that were covered by a slit and 
bright enough for a spectrum to be extracted. In the following, 
we only consider those objects that are either detected in at least 
five filters to attempt SED modelling, or fall in one of the slit 
positions. In total, 11 objects within 100 kpc of the afterglow 
position fulfil these criteria (Fig.[2l Table|4l). Among them, 9 fell 
into one of the slits. 

The two brightest objects (B4 and BIO) within the 100 kpc 
radius are late-type stars. The spectrum of B9 exhibits two emis- 
sion lines at 5322 .1 an d 7149.3 A both detected at ~ 6cr con- 
fidence level (Fig. lA.lb . We identify these lines as [On] A?>121 
and [O III] /15007 at a common redshift of z = 0.428. 

The spectrum of object B2 (Fig. lA.2l) has a continuum break 
at ~ 6750 A and absorption lines at 6660.2 and 6719.5 A. These 
features are consistent with the B aimer break and Can K&H 
absorption at a redshift of z = 0.693. The S/N ratio of the 
VLT/F0RS2 spectra of B3 and B6 are too low for redshift de- 
termination Q Both objects are detected in the nine filters, see 
Fig|A3]and Fig.|A31 The best fits to their SEDs suggest that B3 
and B6 are galaxies at z = 0.67^^^^ and 1.46^^4^ respectively. 
Object B5 was detected in five filters from 473 to 766 nm (Fig. 
IA.4b . The SED is best described by a young and small galaxy at 
z = 0-9^Q 2. We caution that the solution is not unique because 
the galaxy is only detected in five filters. The redshifts of B2, 
B3, B5 and B6 match the redshifts of the strong Mgii absorbers 
at z = 0.6915 and 1.4288 towards GRB 050820A, see TableO 

The continua of objects B7 and B8 are visible in the 2D 
spectrum. Due to seeing the spectrum of B8 is not resolved 
into the two objects B8N and B8S as seen in the HST image. 
The extracted ID spectra of B7 and B8 have a low S/N ratio 
and th us their redshifts cannot be determined. Recently, IChenI 
(l2Q12h observed the field with the IR echellete spectrograph 
FIRE mounted on the Magellan telescope. The slit was oriented 
to cover object B 7 and the compound B8. Based on these ob- 
servations, Chen |2012|) found that the galaxy B7 and the host 
complex B8 are a group of galaxies at z - 2.613, in contrast to 
IChen et al.l (12009) who suggested that B7 and B8S are interact- 
ing galaxies forming the strong Mg ii absorber at z = 0.692. 

Although we elucidated the nature of all objects within 3V7 
that are visible in the HST images, we did not identify a possible 
galaxy counterpart of the intervening sub-DLA. We also did not 
find any candidate brighter than F625W = 26.7 (26.6) mag with 
an maximum impact parameter of 50 (100) kpc. 
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Fig. 5. Fit of the intervening hya absorption features towards 
galaxy D7 (zabs = 3.085; logA^(Hi) = 20.7 ± 0.2) and towards 
GRB 070721B (zabs = 3.0939; logA^(Hi) = 20.1 ± 0.3), and the 
Lya absorption feature in the DLA galaxy D3 (zabs = 3.096; 
logA^(Hl) = 21.3 ± 0.2). The fit of D3 nicely shows the Lya 
emission in the red part of the trough. The afterglow spectrum 
(D4) was acquired with VLT/F0RS2, while the displayed spec- 
tra of D3 and D7 were extracted from the VLT/FORS 1 data. The 
fit is shown with a solid line, while Icr errors are displayed with 
the shaded region. 



3.3. GRB 050908 

The afterglow spectrum of GRB 050908 (z = 3.3467), obtained 
with VLT/FORS 1, revealed an intervening DLA at z = 2.6208 
with log A/^ (H l) = 20.8 ± 0. 1 , as derived from the Voigt profile fit 
shown in Fig. [3] We thus revise the conclusion by Fynb o et al.l 
(^009) that this intervening absorber is a sub-DLA. This DLA 
is peculiar for its strong metal lines (£ Wrest (Si 11/11526) = 
2.24 ± 0.06 A) that are stronger than those of most intervening 
DLAs (Fynbo et al. 2009, their Fig. 11). In addition, the after- 
glow spectrum comprises of an intermediate Mg ii absorber at 



^ The spectrum of B3 shows a prominent emission line at ~ 6310 A. 
The bright [Oi] sky emission line at 6300.3 A partly overlaps with this 
feature, making the identification ambiguous. If this is indeed an emis- 
sion line of galaxy B3, it is likely [O ii] A3121 redshifted to z ~ 0.693, 
coinciding in redshift with a strong Mg ii absorber. 



z = 1.5481. The properties of both absorbers are summarised in 
Tables[2]and[3 

At the redshift of the intervening Mgii absorber and the 
DLA, an impact parameter of 100 kpc corresponds to an angu- 
lar distance of 12! and 12V6, respectively, including the uncer- 
tainty in the afterglow position of 0'.'3. The GRB field and the 
galaxies within this impact parameter are shown in Fig.|4l 

An afterglow spectrum was acquired with VLT/FORS 1 
(PA = 0?0) and with Gemini-N/GMOS (RA = IW.O). In ad- 
dition, a further spectrum was obtained with VLT/FORS 1 (PA = 
53? 8) two years later (Table [TJ. Both FORSl spectra cover sev- 
eral faint objects (R ~ 26 mag) within the 100-kpc radius (Fig.lH 
Table |4]). Apart from the afterglow, we only detect a very weak 
emission feature at ~ 4512 A at the position of CI. Assuming 
that this is a genuine emission line, the line is likely hya red- 
shifted to z ~ 2.71, and not [O ii] because of the non-detection 
of U/3 and [Oiii] ^5007. The Gemini spectrum (RA = 110?0) 
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X = 443.0 nm X = 498.0 nm X = 5^ 6.3 nm X = 533.8 nm X = 546.6 nm X= 563.0 nm X = 570.9 nm 

z- 2.644 z = 3.096 z = 3.096 z = 3.096 z = 3.096 z = 3.631 z = 3.096 



Fig. 7. The VLT/FORSl 2D spectrum obtained after the optical afterglow of GRB 07072 IB faded. The photometric counterparts of 
the different traces are labelled according to Fig. [6] and Tabled The host (D4) is only visible via its Lya emission. The afterglow 
spectrum of GRB 07072 IB would be located at the position of D4. Prominent absorption and emission lines of the galaxies and 
intervening absorption-line systems are highlighted and color coded. The intervening DLA present in the afterglow spectrum is the 
very bright galaxy D3. The spectrum of the Lya emitter D7 displays an intervening DLA as well. 




Fig. 6. Similar to Fig. ffl but for GRB070721B (z = 3.6298). 
Object D4 is the host galaxy. Object D3 is the galaxy counterpart 
to the intervening DLA (z = 3.094) towards GRB 070721B. The 
impact parameter of the DLA galaxy is V\1.9 kpc). The line 
of sight of object D7 also traverses a DLA at z = 3.085. The 
distance between objects D3 and D7 is 21V6, corresponding to 
a projected distance of 167.4 kpc at z = 3.09. The inset is a 5'' 
zoom-in on the position of D3. The afterglow position is marked 
by the cross-hair. 



does not cover any of the faint objects visible in the /?-band im- 
age and we do not detect any emission line in these data. 

Since the field is quite crowded and the amount of spectro- 
scopic and multi-filter data are limited, we can only place a shal- 



low upper limit on the brightness of the galaxy counterpart of 
the both absorption-line systems. Assuming a maximum impact 
parameter of 50 or 100 kpc, the galaxy counterpart cannot be 
brighter than 25 mag and 24.7 mag in the R band. 



3.4. GRB 070721 B 



The VLT/F0RS2 afterglow spectrum of GRB070721B (z = 
3.6298) revealed an intervening absorber at z = 3.0939, a sub- 
DLA with strong metal absorption lines Wrest (Si ii /1 1526) = 
1.71 ± 0.54 A; Table [2]). Figure [5] displays the Lya absorption 
profile; the Voigt profile fit gives log A^(Hl) = 20.1 ± 0.3. In 
Figure [6] we show the field of view and highlight the maximum 
impact parameter of 100 kpc, corresponding to 13V3 at z = 3.09, 
including the uncertainty in the afterglow localisation of 0'.'4. 

Four months after the GRB, an additional spectrum with 
PA = -118? 5 was acquired with VLT/FORSl. The spectrum 
covers several objects within the 100 kpc radius (Table |4]). The 
2D spectrum shown in Fig. [7] re veals that this is a very complex 
and perplexing line-of- sight. The galaxy D7, which was also de- 
tected in the 2D- spectrum of the afterglow, is a bright Lya emit- 
ter at z = 3.615 - very similar to the redshift of the GRB. Object 
D4 is the host galaxy of GRB 070721B and is only detected by 
its Lya emission line in the 2D spectrum (Fig.jT]). The emission 
line redshift (zem = 3.631) diff'ers slightly from the absorption- 
line redshift (zabs = 3. 6298), derived from metal li nes in the af- 
terglow spectrum (see iMilvang- Jensen et al ] l2012l for a detailed 
discussion). In the spectrum of the galaxy D3 we detect several 
strong metal lines as well as a strong Lya absorption line at a 
common redshift of z = 3.096. This galaxy has an impact pa- 
rameter of only 1 ''relative to the GRB position. Hence, D3 must 
be the galaxy counterpart of the DLA seen i n the GRB after- 
glow spe ctrum This was also suggested by Fynbo et al.' (20091) 
and Chen et all (ib09). Intriguingly, D7 also has an intervening 
strong Lya absorption line at a very similar redshift (z = 3.085); 
we derive a H i column density of 20.7 ± 0.2 from the Voigt pro- 
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file fitting, as shown in Fig. [50 However, we do not detect metal 
lines associated with this intervening absorption line sy stem0 At 
that redshift the angular separation between D3 and D7 of 2(7/7 
translates into a transverse distance of 161 kpc. The weak trace 
below D7 is from the galaxy D6. This spectrum displays a sin- 
gle strong emission line. It is very likely Lya at Zem = 2.642 and 
not [011/13727], because of the lack of UJ3 and [Oiii] /15007 
in emission. The spectrum of D2 shows that this is also a Lya 
emitter at a very similar redshift to that of D6, Zem - 2.644. 
In the spectrum of the afterglow we detect absorption lines from 
Lya, C II /II 334 and the C iv doublet (unresolved) at Zabs = 2.655 
albeit a low S/N ratio. The velocity distance between the absorp- 
tion and emission line redshifts is ~ 900 km s"^ 

The region around Lya in the spectrum of D3 seems to be a 
superposition of an emission line on top of a broad absorption 
line (see Figs. [51 0. Such a feature has be e n obse rved in several 
high-z galaxies (e.g. 'Pettini et al.lll998allbl l2000l) . We estimate 
an H I column density of log (H l) = 21 .3 ± 0.2. 

4. Discussion 

4.1. Detected galaxy counterparts 

4.1 .1 . DLA galaxy towards GRB 070721 B 

In Sect. 13.41 we showed that the galaxy counterpart of the inter- 
vening DLA (object D3, zd3 = 3.096 ± 0.003) has an impact 
parameter of 7.9 kpc and an extinction-corrected /?-band magni- 
tude of 24.41 ± 0.02 mag (Fig.[6l Table[4]). The co-ordinates are 
RA(mOO) = 02^2^33^.018 and Z)£C(J2000) = -02°ir4a'99 
with an uncertainty of 0'.'4 in each coordinate. We denote the 
object as DLA J0212-0211 in the following. Up to now, only 
o ne further sub-DLA/DLA galaxy was identified beyond z > 
3. iDjorgovski et al.l (Il996l) found the galaxy counterpart (R = 
24.8 + 2 mag ) to the z = 3.15 sub-DLA (logA^(Hi) = 20.0; 
iLu et al.|[T993l) towards QSO B2233-hl31. Compared to typical 
z - 3 galaxies, both objects are comparable to the brightness of 
an galaxy (R^ = 24 AS ± 0.15 mag; S teidel et al.|[T9991 

At low-z (z < 1), the overwhel ming majority (75 %) of 
sub-DLA/DLA galaxies summarised in Peroux et al. (2011) are 
fainter than both z> 3 sub-DLA/DLA galaxies. The mean lumi- 
nosity of that sample is 0.66 L* and the median is 0.32 L.,. The 
sample is not tightly distributed around the mean and the bright- 
est galaxy in that ensemble reaches 2.8 L^,. A proper comparison 
with high-z (z > 2) DLA galaxies is not possible. To date, the lu- 



^ There are some uncertainties in the Lya profile fits for D3, D4 and 
D7. Because of the resolution and the quality of the spectra, the damped 
wings are barely visible (see Fig. [5]). We cannot rule out the possibility 
that a blend of narrow Lya lines mimics the strong absorption feature. 
It is not very likely, because it would require that D3 and D7 have a 
comparable density of H i clumps at almost identical redshifts. Strictly 
speaking, the stated column densities are only upper limits. The fit does 
not rule out an extremely large Doppler parameter (b » 100 km s~^), 
implying a reduction of the H i content by a factor of 100. This remark 
has in particular to be kept in mind for the strong Lya absorption feature 
of the intervening absorber towards D7. 

^ We obtain the following 3<t upper limits from the VLT/FORSl 
spectrum, assuming an aperture of 1665 km s~^ (twice the FWHM 
of a Gaussian with a Doppler parameter of b = 500kms~^): 
£Wrest(SiiU1260) < 3.6 A, ^W^est (0 1/11302) < 2.8 A, 
£Wrest(C 11/11334) < 2.8 A, £Wrest(SiiV/11339) < 2.7 A, 
Wrest (C IV /11548) < 2.7 A, £Wrest(Fe 11/11608) < 3.0 A, and 
EWrcst (Al 11/11670) < 3.3 A. These limits, are on average less stringent 
than the measurements of the detected metal lines in the DLA galaxy 
D3 (see Table[6]). 



minosity of only 4 of 10 sub-DLA/DLA galaxies is known. The 
reason for this is that most of them were detected by their Lya 
emission and have small impact parameters so that they are out- 
shone by the glare o f the quasar . This hampers the determination 
of their luminosity. IChen et al l (l2QQ9 ) concludes that the DLA 
galaxy in the field of GRB 07072 IB is the most luminous DLA 
galaxy at z > 2. The morphology of the galaxy is undisturbed, 
based on the shape in the HST image, inset in Fig. [6] We measure 
an ellipticity of 0.15 and a half-light radius of 1.2 kpc, obtained 
with SExtractor. The non-detection o f a galaxy interaction or 
merger is not surprising. lOverzier et al. I {2010) argued that a de- 
tection is hampered at z ~ 3 because of the reduced physical 
resolution and sensitivity in addition to the general difi&culties of 
observing interacting galaxies. 

The Lya absorption line shows an excess of flux in the 
red part of the absorption profile (Fig. [5} Sect. 13.41) . The 
extracted emission line profile, shown in Fig. [8] is slightly 
asymmetric and peaks at ~ 1000 km s"^ with respect to 
the systemic redshift. We measure a line flux density of 
(2.34 ± 0.25) X 10~|^ ergcm"^ s~^ A ~\ i.e. a line significance of 
9.2 cr. iPettini et al.l (Il998allbl l2000l) reported the detection of the 
same feature with peak recessional velocities b etween -\-400 and 
1100 km s"^ in several z ^ 3 galaxies (see also lAdelberger et aP 
I2003h . They argue that this feature is similar to a P-Cygni profile, 
indicating the prese nce of a gala ctic outflow (also found in star- 
burs t galaxies; e. g. iKunth et al.l J9 98: Gonz alez Delgado et al.l 
[1998 : but see also [Verhamme et al.[ r2006: Lau rsen et al.[[2009[) . 

To further explore the properties of the galaxy counterpart, 
we fit the SED with Le PHARE (Sect. [23]), as displayed in 
Fig. [8] Leaving all model parameters free, except for the red- 
shift, which was fixed to z = 3.096, the SED is best described by 
a young (0.4 Gyr) and close to dust-free galaxy (Ay = 0.3 mag, 
Calzetti reddening; see Table [5] for all best fit values). The SFR 
of 37 Mr?) yr~^ is similar to the that of the z = 3.15 DLA galaxy 
(Dio rgovski et al.[ [19961: [Christensen et al.[[2004bl: [Peroux et alJ 
201 I j) and similar to the SFR of typical t/-band dropout LBGs 
(Gia vafisco[[2002h . 

Because of the small impact parameter of the DLA coun- 
terpart (V), the GRB line-of-sight can be extingu ished by the 
dust in the DLA gala xy if it cont a ins d ust (see iFynbo et al.[ 
2011 for an example). [Fynbo et aP ([2009i) derive the spectral 
slope betwee n the optical and X-ra y bands, y^ox, to be 0.72. 
According to Jakobsso n et al.[ ([20041) . a dust-extinguished opti- 
cal afterglow woul d result in jSg^ < .5, which is not the case 
for GRB 070721B. [Zafar et al.[ ([201 ih obtained a visual extinc- 
tion of Ay = 0.20 ± 0.02 mag, assuming that all extinction is 
attributed to an absorber in the host galaxy of GRB 07072 IB 
This value is typical for GRB host galaxies (Kan n et aD[2010[: 
[Schadv et al.[[2010[:[Zafar et al.[[201 ih . Motivated by the fact that 
some GRBs show negligible host reddening, the presence of a 
foreground absorber cannot be excluded. Assuming the DLA to 
be responsible for the obs erved reddening , we expect that Ay in- 
creases only slightly (Grei ner et al.[[201 U their Fig. 3). Hence, 
we can in this case exclude the presence of significant amounts 
of dust, i.e. Ay < 0.20, in the outskirts of the DLA galaxy. 

In addition, we measure the EWs of absorption lines in the 
DLA galaxy spectrum D3 and of absorption lines imprinted from 
the DLA galaxy on the GRB afterglow spectrum. The EWs and 
column density estimates are listed in Table [6] By comparing 
the two sets of data it is evident that the DLA galaxy shows a 
higher N (H l) than in the afterglow spectrum (see Fig. [5]). This 
trend is also confirmed by all the metal EWs being higher for 
the DLA galaxy, despite the individual large errors. These re- 
sults are easily explained if most of the DLA galaxy light arises 
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Table 6. DLA J02 12-021 1 absorption lines in the DLA galaxy and the GRB 07072 IB afterglow spectrum 



Ion 



Transition 
(A) 



(A) 



\ogN 



[X/H] 



(A) 



logiV 



[X/H] 







Spectrum of the DLA galaxy D3 


D3 in the afterglow spectrum 


Hi 


1215 




21.3 + 0.2^ 






20.1+0.3^ 




Sin 


1260 


3.00 + 0.98 


> 14.32^^ 


> -2.49 


< 2.44 






Oi/Siii 


1302/1304 


3.94 + 1.15 






3.06 + 0.83 






Cii 


1334 


3.97 + 1.11 


> 15.29 


> -2.44 


< 3.57 






Si IV 


1393 


3.41 + 1.01 


> 14.57 




<2.59 






Sin 


1526 


< 5.05 






2.35 + 0.69 


> 14.95^ 


> -0.66 


Civ 


1548/1549 


< 6.13 






4.34 + 0.73 


> 15.03 




Fe II 


1608 


<4.98 






3.02 + 0.62 


> 15.39 


> -0.21 


Mil 


1670 


<7.33 






2.33 + 0.74 


> 13.70 


> -0.85 



Notes. The last column denotes the metal abundance with respect to solar. The restframe EW (Icr errors) and the 3cr detection limits are 
measured assuming an aperture of 1665 km s~\ i.e. twice the full width at half maximum of a Gaussian with Doppler parameter b = 500 km s~^ 
(b) Y)QYi\ed from the Voigt profile fit of the hya absorption, see Fig.|5] The column densities, derived from the EW in the optically thin limit, are 
considered as lower limits because the lines are saturated. 



Table 5. Properties of the galaxy counterpart candidates of the 
strong Mgii absorber towards GRB 05 08 20 A (zMgn = 0.692 
and 1.430) and the DLA galaxy towards GRB 070721B (zdla = 
3.094) 





B2 


GRB 050820A 
B3 B5 


B6 


GRB 07072 IB 
D3 


z 


0.693 


-0.06 


9+01 


1 46+007 


3.096 




87.5 


36.8 


26.4 


42.0 


7.9 


Mb (mag) 


-20.2 


-19.4 


-17.7 


-21.4 


-22.3 


LIU 


0.3 


0.2 


0.03 


0.9 


~ 1 


Dust model 






Calzetti 






Av (mag) 


0.3 


1.6 


0.3 


0.3 


0.3 


Age (Gyr) 


2.0 ±0.2 


1 3+2-3 




5+0-5 


4+0.5 


log M/Mq 


10.3 ±0.1 


Q Q+0.2 


8.4 ±0.2 


10.1 ±0.1 


10 1+0-3 


SFR (Moyr-i) 




7.5 ±2.7 


5+1-2 


1 5+1-5 


07+117 
-"-23 




16.3/9 


7.3/9 


1.5/5 


1.6/9 


5.0/4 



Notes. The table summaries the most important properties of the SED 
fits of the most likely candidates. The redshifts of B2 and D3 were fixed 
to the spectroscopically measured redshift in the SED fit. The impact 
parameters, 6, are calculated from Table ID assuming either z = 0.6915, 
1.430, or 3.096. Th e lun iinosities were com puted using the results by 
iDahlen et all (l2005h and ISteidel etaP ([1999) on EEs; the knee of the 
luminosity functions is at Mb,* = -21.43 Mb - -21.60 and mr - 
24 AS mag at z = 0.692, 1.4288 and 3.096, respectively. We note that the 
stellar mass and its uncertainty wa s de rived assuming a Chabrier IME 
dChabrier et al.l200 Q. see also Sect.Ofor more details on the SED fit), 
"n.o.f." stands for number of filters. 



from the central, i.e. denser, region of the galaxy, while the GRB 
afterglow radiation shone through a more peripheral region of 
the host. 

Finally, the most constraining metallicity estimate for the D3 
galaxy is [Fe/H] > -0.21 (measured from the GRB afterglow 
spectrum). Iron being a very refractory element, a significant 
amount of Fe could be depleted into dust grains and not observed 
in the gas-phase. If this is the case, then the intrinsic metallicity 
would increase to solar or super-solar values, while GRB DLAs 
typica lly show solar or sub-solar metallicities (see iFynbo et aP 
l2009l and references therein). A more likely scenario is a lim- 
ited amount of dust in the DLA galaxy, also supported by the low 
reddening along the line-of- sight (Ay < 0.20) estimated from the 
SED. Thus, the metallicity along the GRB line-of- sight, i.e. 7.9 
kpc away from the galaxy centre, is close to ^ 0.67©. A higher 
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Fig. 8. Spectrum (middle) and the SED (bottom) of the DLA 
galaxy D3 in the field of GRB 070721B. In the middle panel, 
the error spectrum is overplotted and several absorption lines are 
marked. Regions of strong telluric features (atmosphere trans- 
parency < 20%) were not used in the spectral analysis and are 
marked by a small box at the top. The top panel shows zoom-ins 
on the individual absorption lines. In the bottom panel, the ob- 
served extinction-corrected data points are shown as circles with 
error bars. The curve represents the best fit to the observed SED. 
The model predicted magnitudes (squares) are superposed. The 
gray area highlights the interval that is covered by the spectrum 
above. 



metallicity could possibly be found in the central regions of the 
DLA galaxy, given the likely metallicity gradient, but cannot be 
constrained from our dataset due to absorption-line saturation. 

In conclusion, the properties of the DLA galaxy D3 are 
very typical for LBGs, except that the halo of the galaxy har- 
bours a large amount of neutral hydrogen. This is consistent 
with the model that DLAs are gaseous halos of faint high-z 
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LBGs (jFvnbo et al ] |1999I: lM(z)ller et al.ll2QQ2l: iFvnbo et al.ll2008l: 
iRafelski et al.ll2011h . 



Table 7. Limiting magnitudes of the galaxy counterparts of the 
intervening absorption line systems 



4.1 .2. Mg II absorbers towards GRB 050820A 

In addition to the intervening sub-DLA towards GRB 050820A, 
the afterglow traverses two strong Mgii absorbers at z = 0.6915 
and 1.4288 (see Table [3]). We find four galaxies within 100 kpc 
that are consistent with either z = 0.6915, namely B2, B3 and B5 
(Figs. |A21|A31 (All Tabled]), or z = 1.4288, namely B6 (Fig. 
IA.5I Table OH The properties of the best-fit galaxy templates 
are summarised in Table [33 

The impact parameters vary between 26 and 88 kpc for the 
candidates of the Mgii absorber at z = 0.6915. Galaxy counter- 
parts with an absorption cross-section of more than 88 kpc are 
rare; for instance ChenetaL (12010.) used a set of 94 galaxies 
((z) = 0.24) that are located up to 120 kpc around 70 back- 
ground quasars (zqso > 0.6) to study which kind of Mg ii ab- 
sorber could be found in the quasar spectrum arising from these 
galaxies. They did not find any strong Mg ii absorber in a quasar 
spectrum that is related to a galaxy at a distance of more than 
~ 40 kpc from the line of sight of the quasar. Hence, the early- 
type galaxy B2 is probably not the galaxy counterpart of the 
absorption-line system. The impact parameters of B3 and B5 are 
typical for strong Mgii absorbers (iChen et al.l 12010), though at 
the upper end of the impact parameter distribution. The proper- 
ties of B3 and B5 are very diff'erent from each other (TableO, B3 
is a massive galaxy with a high S FR, while B5 is a young low- 
mass galaxy with low S FR. These diff^erences are also reflected 
in their luminosities; adopting the nDahlen et al.l (l2005h LF, B3 
is 0.2 L*, and B5 is 6 times fainter (0.03 L^). Both values are in 
the range of the observed luminosities of the IChen et aL 
sample, but in their fainter half. 

Galaxy B6 is the most likely galaxy counterpart of the strong 
Mgii absorber at z = 1.4288. The impact parameter of 42 kpc 
(4V9; Table O is similar to the impact parameter of the galaxy 
counterpart candidates of the Mg ii absorber at z = 0.6915. The 
properties of the best-fit galaxy template, displayed in Table 
[51 show that it is a fa irly bright (0.9 L^,) galaxy with S FR = 
1.5Moyr-i (Fig.[A3]). 

Without additional spectroscopic data it is not possible to 
decide which of two candidates is the galaxy counterpart of the 
strong Mgii absorber z = 0.6915 and if B6 is the galaxy coun- 
terpart to the Mgii absorber at z = 1.4288. Chen' ('2012') places 
an upper limit on the brightness of the galaxy counterpart of 
F115W = 27.5 mag, i.e. fainter than a 0.03 galaxy, if the im- 
pact parameter is less than 3V5. 

4.2. Limits on galaxy counterparts 

In the previous sections we discussed the properties of the de- 
tected DLA galaxy towards GRB 07072 IB and of the galaxy 
counterpart candidates of the two strong Mg ii absorbers towards 
GRB 050820A. In this section we present the limits, such as lu- 
minosity and S FR, on the non-detected galaxy counterparts. The 
limits on the luminosities s hown in Table[7]are re lated t o the knee 
of the lDahlen et al.l (l2005l) and the .Reddv etaP (12008) LFs. The 



The photometric redshift of B6 also matches that of a weak Mg ii 
absorber at z = 1.6204 towards GRB 050820A. It is unlikely that the 
galaxy is the counterpart of the weak Mg ii absorber, because the corre- 
lation length of weak Mg ii absorbers is ~ 2 kpc (E llison et al . 2004). 

^ ^ We remind the reader that the solution of B 5 is not unique, because 
it is only detected in five filters. 



Absorber Zabs 


FllSWuL (mag) for 
6* < 50 (100) kpc 


L/L* 


Fynbo 
model 


GRB 050730 (zgrb = 


3.969) 






DLA 3.56439 
sub-DLA 3.02209 
Mgii 2.25313 
Mgii 1.7743 


> 25.7/ > 24.6 

> 25.7/ > 24.6 

> 26.5/ > 25.1 

> 26.5/ > 25.1 


<0.5/ < 1.5 

< 0.4/ < 1.1 

< 0.2/ < 0.5 

< 0.1/ < 0.5 


< 0.04 

< 0.04 


GRB050820A (zgrb 


= 2.615) 






sub-DLA 2.3598 


> 26.7/ > 26.7^ 


< 0.1/ < 0.1 


<0.1 


GRB 050908 (zgrb = 


3.3467) 






sub-DLA 2.6208 
Mgii 1.5481 


> 25.0/ > 24.7^ 


< 0.4/ < 0.6 

< 0.3/ < 0.4 


>0.5 



Notes. The lim iting magnitudes are cor rected for Galactic extinction. 
We adopted the lDahlen et"an 12005) and' Reddv et~ (2008) LFs to re- 
late an absolute magnitude with the knee of a corresponding LF, L*: 
1.59 < z < 1.9 : M, = -20.69 + 0.26 mag (U-band), 1.9 < z < 
2.7 : = -21.01 + 0.38 mag (1700A), 2.7 < z < 3.4 : = 

-20.84 + 0.12 mag (l700 A). We extended the range of the LFs from 
z = 1.59 to 1.50 and from z = 3.4 to 3.56, because the change in the 
is smaller than its uncertainty. To compute the absolute magnitude we 
used: M - m - DM(z) + 2.5 log (1 + z), where DM (z) is the distance 
modulus. Different filters than F775W were used if a filter was closer to 
the central wavelength of the LF than the F775W band. In the column 
"Fynbo model" we display the limits on the luminosity of DLA galaxies 
based on the model by Fynbo et al. (2008). The F625W magnitude 
was used. The R-bmd magnitude was used. 



galaxy counterparts are fainter than a 0.5 and 1.5 galaxy in 
all fields, if the impact parameter is 50 and 100 kpc, respectively. 
The most stringent upper limit of 0.1 can be placed on the 
galaxy counterparts of th e intervening sub-DLA towards GRB 
050820A. For this field, IChenI (l2012l) reports that the galaxy 
counterpart is fainter than F775W = 27.5 mag, corresponding 
to 0.06 L^, if the i mpac t parameter is smaller than 3V5. 

iFynbo et al.l (IT999I) suggested tha t DL As are the g aseous 
halos of LBG (see also iM^ller et al.l [20021: [Fvnbo et al.[ [2008[: 
[Rafelski et al.[ [201 H) . To confront our measurements with this 
model, we use the predictions on the brightness of galaxy coun- 
terparts by [Fynbo'e t al. (2008). In their model, simple and con- 
strained scaling relations of galaxies in the local Universe are 
used to predict the luminosity of intervening DLAs based on 
their metallicity and their impact parameter. These relations are 
strictly speaking only valid for z ~ 3 DLAs. We also adopt this 
model for DLAs at lower redshift and for sub-DLAs. The main 
diflTerence between a sub-DLA and a DLA is a larger covering 
fraction, i.e. impact parameter, if DLAs and sub-DLA are drawn 
from the same population of galaxies. Based on the measured 
metallicities, we estimate the brightness and the impact parame- 
ter of the unidentified DLA galaxies towards GRBs 050730 and 
050820A. For G RB 050908, we used ^ Wrest (Si 11.11526) as a 
metallicity proxy (Prochask a et al.[[2008[) . 

The estimated brightnesses, shown in Table [T] indicate that 
most galaxy counterparts evaded detection because of their in- 
trinsic faintness. This result is not surprising. The majority of all 
z ~ 3 DLA galaxies are expected to be fainter than ~ 29. 5 mag in 
the 7?-band, i.e. L < 0.012 L„ in the Fynbo etalJ ([2008 ) model. 
The model clearly rules out large impact parameters. Even the 
larger covering fraction of sub-DLAs does not imply impact pa- 
rameters of more than 2-3'' (15.7-23.5 kpc at z ~ 3). The in- 
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Table 8. Absorption and emission lines of QSO J1408-0346 and 
the strong Mg ii absorber at z = 1.77425 towards it. 



^obs 

(A) 


1 

^rest 

(A) 


opcuics 




^ »^»^rest 

(A) 


3658 


911.20 


Ly-limit 


3.01 




^ 4168.8 


^ 1033 


Lyyg + Ovi 


^3.04 




4888 


1215.80 


Lyof 


3.02 




4992.72 


1239.42 


Nv 


3.028 




5621.30 


2026.14 


Zn II 


1.77425 


< 0.73 


6208.75 


1548.20 


Civ 


3.010 




6210.62 


1548.20 


Civ 


3.011 




6219.37 


1550.77 


Civ 


3.010 




6220.87 


1550.77 


Civ 


3.011 




6229.35 


1545.86 


Civ 


3.030 




6503.52 


2344.21 


Fell 


1.77425 


<0.49 


6610.51 


2382.77 


Fe II 


1.77425 


0.84 + 0.15 


7145.59 


2586.65 


Fell 


1.77425 


<0.54 


7213.61 


2600.17 


Fe II 


1.77425 


0.62 + 0.19 


^ 7674.8 


^ 1910 


Sin, Ciii/Feii 


3.018 




7756.18 


2796.35 


Mg II 


1.77425 


1.24 + 0.34^ 


7777.52 


2803.53 


Mg II 


1.77425 


1.25 + 0.19 



Notes. The table summarises the most important absorption and emis- 
sion lines of the quasar and the strong Mgii absorber at z = 1.77425. 
We show only those absorption lines that are redward of the quasar Lyof 
emission line. Emission lines are displayed in bold, 
(a) -pj^g restframe EW^ (Icr errors) and the 3cr detection limits are mea- 
sured assuming an aperture of 333 kms~\ i.e. 2 x FWHM a Gaussian 
with Doppler parameter b = lOOkms"^ The EW of the lines of the 
Mg II doublet are similar, indicating that the lines are saturated. 



tervening DLA towards GRB 050908 possibly has a ~ 0.5 
galaxy counterpart. 

We find several objects in the field of GRBs 050730 and 
050908 that fulfil both criteria, being very faint and having a 
very small impact parameter (Table |4l). There are several objects 
with a brightness of F115W ~ 28.4 mag within 2! 5 of the af- 
terglow position of GRB 050730. If one of these objects is the 
galaxy counterpart to the sub-DLA, it is a 0.03 L* galaxy. In the 
field of GRB 050908 we find four objects ranging in brightness 
from R = 25.3 to 27 mag within 4'' of the afterglow position, 
corresponding to L < 0.4 L.,. To elucidate the nature of these 
objects, more data is required. 

4.3. Quasar radiation field 

In Sect. 13.11 we alluded to the serendipitous coincidence of the 
sub-DLA (z = 3.02) towards GRB 050730 and the QSO that is 
17! 5 south of the GRB position. Before we discuss their con- 
nection, we present general properties of the quasar. 

After subtracting the bright R = 18.2 mag foreground star 
(Sect. l2.2l Tablelll), we measure an 7?-band magnitude of 20.78 ± 
0.05 mag (corrected for Galactic extinction. Table |4]). The co- 
ordinates of the quasar are 7?A(J2000) = 14^08^17^35 and 
DECimOO) = -03°46'35r35, with an uncertainty of (y!4. 
Hereafter the QSO is denoted as QSO J1408-0346. Figure [9] 
shows the quasar spectrum secured with X- shooter on 18 April, 
2010. Table [8] summarises the detected quasar emission and ab- 
sorption lines; the mean emission line redshift is Zem = 3.023 ± 
0.007. The emission lines are substantially narrower i n com- 
parison to the quasar composite spectrum by lParis eTal ](|2011), 
which is overplotted in red. In addition to the emission lines, we 
detect Civ and Lya in absorption at Zqso (Table [U). The narrow 



C IV absorption lines consist of two blueshifted velocity compo- 
nents at 1420 and 1497 km s"^ with respect to the maximum of 
the C IV emission line. 

The proximity of the GRB l ine-of - sight can constrain the 
quasar radiation field. iFox et al.l (l2QQ7l) reports that DLAs and 
sub-DLAs bearing Nv are extremely rare, and that several of 
these have a nearby QSO coincident (> 33%) at the ( sub)- 
DLA redshift (see also Fox et al. 20071: lEUison et"aDl201Ql) . We 
do not find evidence for Nv in absorption in the GRB after- 
glow spectrum. Besides this, the ionising UV radiation field of 
a bright quasar can also lead to visible underdensity in the Lya 
forest as the redshift approaches the qua sar emission redshift . 
This is referred to as the proxim ity eff'ect (ICar swell et al.ll 19821 
lTvtlejll987l: iBaitlik et al.lll988l) . In addition, if the quasar radi- 
ates isotropically, the line of sight of GRB 050730 might show 
an underdensity in the hya forest around the redshift of the fore- 
ground quasar (i.e. the transverse proximity eff'ect; Bajtlik et all 
1988). At the projected distance of 137 kpc, the QSO with an 
inferred Lyman-limit flux density of 9 ju Jy exceeds the Hi UV 
backgrou nd at z ~ 3 by a factor of ~ 170 ( Guimaraes et al. ^20071; 
Dan'Agli oetani2009h . In this case, the hya forest at the redshift 
of the QSO is dominated by the sub-DLA so it impossible to in- 
fer if there is an efl'ect of the QSO on the IGM along the GRB 
afterglow line-of- sight. 

4.4. Correlated structures 

The spectra of GRB 050730 and the QSO 17^5 south of it (Fig. 
[TJ share a common wavelength interval, as do GRB 07072 IB 
and the bright galaxy D7 20^.7 SW of it (Fig. O. These allow 
us to identify correlated structures in both fields. We are in par- 
ticularly interested in absorbers with metal lines, because they 
are related to galaxies. Their detections provide information on 
galaxy impact parameters and their clustering. 

Correlated structures have been found in several quasar 
pairs on a l ength scale of a few k pc up to ~ 100 kpc (e.g . 
Smette etal.lll99l iPetrv et al.llT998l: lEUison et al.ir2004t l2007h . 
and on a length sca le of several hundreds of k pc up to a few Mpc 
("Francis &Hewett'"1993"; 'Francis et al."1996':' Fvnbo et al.ll2QQ3l: 
D'Odorico et al. 2002; Coppolani et al. 2006). In the literature, 
two diff'erent scenarios are discussed to explain correlated struc- 
tures on these diff'erent length scales: a) the lines of sight probe 
the halo of one galaxy, or of a galaxy and a satellite galaxy, and 
b) the diff'erent lines of sight probe diff'erent galaxies that belong 
to a group of galaxy and are part of a large scale structure. Here 
we argue that the latter interpretation of correlations being due 
to large scale structure explains better our observations of the 
absorbers in both GRB fields. 

4.4.1 . The field of GRB 050730 

The lines of sight of GRB 050730 and QSO J1408-0346 are 
separated by 17V 5. We identify a Mgii at z ~ 1.774 in both of 
them, as well as the sub-DLA towards the GRB and QSO J1408- 
0346 that are at almost identical redshifts (Zsub-DLA = 3.02209, 

^QSO,em = 3.023). 

The coincidence of the sub-DLA and the QSO is in line with 
evidence for anisotropic emission of QSOs, which seem to have 
cleared their surroundings prefer entially towards us rather tha n 
perpendicular to the line-of- sight (iHennawi & P rochaskal l2QQ7h . 
At le ast it suggests that QSOs reside in over-de nse regions 
(Rol linde et al.ll2QQ5l: iFaucher-Giguere et al.ll2QQ8h . The angu- 
lar distance of 17V 5 between the sub-DLA and the QSO trans- 
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Fig. 9. X-shooter spectrum of QSO J1408-0346 (bottom). The spectrum beyond 1.8 yum is not shown due to uncertainty in the flux 
calibration. The observed and rest-frame wavelengths are shown below and above each panel, respectively. The error spectrum is 
shown in grey. Regions of strong telluric features (atmosphere transparency < 20%) were not used in the spectral analysis and 
are marked by small boxes a t the top of each panel (NIR: Gemini observatory). The red line is the average z = 3 QSOs spectrum 
taken from Paris et al.l (l201ll) . Prominent absorption and emission lines from the quasar are indicated by vertical solid lines, and the 
absorption lines from the strong Mgii absorber at z = 1.77425 by vertical dashed lines (see Table [8]). Zoom-ins on the absorption 
Iji^es of the Mg ii absorber are shown in the top panel. They were normalised with the QSO composite spectrum. 
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lates into a transverse distance of 137 kpc at z - 3.022, which is 
several times larger than the typical sub-DLA impact parameter 
of ~ 39 k pc in Peroux et al. (2011) and Krogager et al. (2012). 
lEUison et al.l (l2010l) argue that proximate DLAs, i.e. DLAs 
within a comoving distance of 42 Mpc (Av = 3000 km s"^) 
around QSOs, are not associated with the Q SQ host, but rathe r 
sample overdensities around it (see also Russ ell et"an 12006). 
Therefore, the intervening sub-DLA towards GRB 050730 could 
be an overdensity in the vicinity of the QSO rather than a part of 
the massive halo of the quasar host galaxy. 

The common Mg ii absorbers towards GRB 050730 and the 
0SO_^re_at ZQRB.Mgii = 1.7731 and ZQso,Mgii = 1.77425. 
IVergani et al.l (|2Q09>) reported a restframe Mg ii equivalent width 
of 0.93 ± 0.03 A and lProchaska et al.l (l2007ah reported that the 
doublet is saturated so that the absorber can be placed in the 
strong Mgii absorber category. For the intervening Mgii ab- 
sorber towards the QSO, we measure £^Wrest(2796 A) = 1.24 ± 
0.34 (saturated; Table[8l). In addition, we detect Znii A 2026 and 
several Fe ii absorption lines of that absorber at the same redshift, 
summarised in Table [S] At the redshift of 1.774, the angular dis- 
tance of 17! 5 translates into a projected distance of 150 kpc, 
while their redshift implies a velocity distance of 124 km s"^. 
IChen et aD (l2010l) reports a typical impact parameter of around 
~ 30 kpc for Mgii absorbers (see also Smette et al. 1995), ar- 
guing against the hypothesis that both strong Mg ii absorbers are 
associated with the same galaxyQ Finally, we note that there is a 
high-column density absorber at z ~ 2.98 in the afterglow spec- 
trum, which has no counterpart in the spectrum of the quasar. 

4.4.2. The field of GRB 070721 B 

In Sect. 13.41 we identified object D3 as the galaxy counter- 
part of the intervening DLA in the afterglow spectrum of GRB 
07072 IB. In addition, we identified another intervening DLA to- 
wards the Lya emitter D7 that is at a similar redshift and has a 
similar Hi column density to the DLA galaxy, but is 160.7 kpc 
from the location of D3 (Figs. [51 0. It is unlikely that the inter- 
vening DLA towards D7 is gravitationally bound to the DLA 
galaxy D3 because of this large impact paramet er. The aver- 
age impact parameter of DLAs is 13 kpc, although [Francis et"aD 
(1996) also report a candidate a counterpart with impact param- 
eter of 182 kpc. 

D3 and the intervening absorber towards D7 appear more 
similar to the pairs of intervening absorption-line systems dis- 
cussed in Ellison et al. (2007)0 They found found a pair of sub- 
DLAs at z = 2.94 and a pair of DLAs at z = 2.66 towards a 
binary quasar, where the individual absorbers of each pair are 
separated by ~ 100 kpc. Based on simulations, they showed that 
the presence of a large scale structure is more likely than the 
diff'erent lines of sight probing the halo of two massive galaxies 
at z = 2.66 and 2.94, respectively. In this case, the probabil- 
ity increases from p < 10"^ to 0.01. We find several correlated 
structures on a length scale of ~ 160 kpc in the fields of GRBs 



The galaxy counterpart of the Mg ii absorber towards the QSO is 
fainter than 26.5 (25.1) mag in F115W (corrected for Galactic extinc- 
tion) for an impact parameter of 50 kpc (100 kpc), assuming that the 
galaxy counterpart is neither in the glare of the QSO, nor a t the posi- 
tion of the foreground star. Adopting the Dahle n et"aLl (I2005D LP, the 
limiting magnitudes correspond to L < 0.04 L* and L < 0.16L*. 

The galaxy counterpart of the intervening DLA towards D7 is 
fainter than 25.2 mag in F115W, assuming an impact par ameter of 50 
to 100 kpc, corresponding to L < 0.7 L*, using the LP in lReddv et al.l 
(i2008i) . 



050730 and 07072 IB. These are most likely further examples 
for a group environment of intervening absorption-line systems. 

5. Conclusions 

The aim of our work is to detect galaxy counterparts of high-z 
intervening sub-DLAs and DLAs. In contrast to previous stud- 
ies, we use GRBs to have a clear view on the region that is usu- 
ally outshone by the glare of a quasar. Since the launch of the 
Swift satellite, seven intervening sub-DLAs and DLAs have been 
found towards six GRBs. Among them four lines-of- sight have 
sufficient photometric and spectroscopic data to study very faint 
objects. 

In our study we successfully detected the DLA galaxy, caus- 
ing the intervening DLA absorption towards GRB 070702 IB, V 
fror n the afterglo w position, as suggested by Chen et al. ( 2()0i) 
and lFynbo et al.l (|2009). However, these authors did not present 
any spectroscopic evidence for their inference. In fact, the DLA 
galaxy would have been almost impossible to detect by direct 
imaging if the background source would have b een a 19-mag 
quasar. This underlines the argument by Jakob ssonetaD (12004)) 
that studies of the galaxy counterparts of intervening absorption 
line systems towards QSOs can be aff'ected by misidentifications. 
Hence, proximity is not sufficient for an association. For instance 
in the case of GRB 050820A, the extensive photometric and 
spectroscopic campaigns allow us to successfully rule out all ob- 
jects brighter 26.2 mag in F625 W-band within 3V7, correspond- 
ing to L < 0.1 L:, at z = 2.3598. Assuming that su b-DLAs and 
DLAs are LBGs weighted by their H i cross-section. lFynbo et al.l 
(2008) showed that the overwhelming majority of intervening 
DLA galaxies are expected to be fainter than this, naturally ex- 
plaining the non-detections. Even with the largest telescopes, it 
is difficult to detect and elucidate the nature of such faint objects. 
On the other hand, there are cases where intervening absorption- 
line systems w ere successfully associat ed with galaxies at larger 
distances (e.g. lFrancis & Hewettll 19931) . however they are just a 
minority. 

The disadvantage of our approach is the vast amount of ob- 
serving time required. Deep multi-filter observations exist for 
several GRBs. However, similar extensive spectroscopic cam- 
paigns do not exist for other G RB fields. In tho se cases, one 
relies on SED fitting techniques. iRao et al.l (1201 ih showed that 
this can indeed be successful, in particular for very faint objects 
where only emissi on-line spectroscopy is feasible. Recently, 
iPeroux etaP (l201ll) performed a survey for galaxy counterparts 
of intervenin g DLAs using the integ ral field unit VLT/ SINFONI 
(see also Chri stensen et al.i r2004a: B ouche et al.ll2012h . This ap- 
proach could be complementary to the classical strategy of per- 
forming very deep imaging campaigns and spectroscopic follow- 
ups on candidates, or the use of narrow band and broad band 
strategy to search for emission line objects (e.g. lVreeswijk et all 
2003). 

The DLA galaxy (z = 3.096 ± 0.003) in the field of 
GRB 07072 IB is the most luminous high-z DLA galaxy known 
dChen et al.l2009h . and is almost as distant as the highest-redshift 
DLA galaxy known so far at z = 3 . 1 5 . The number of high-z sub- 
DLA and DLA galaxies increases thus from nine to 10. DLA 
J0212-0211 is very metal rich, [Fe/H] > -0.21. It does not dif- 
fer from normal LBGs, with the exception of the large amount 
of H I in its halo. This supp orts a model in which DLAs can be 
gaseous halos of LBGs (Fvnbo et al.Jll999|; lM0ller et al.l 120021: 
.Fvnbo et al. 2008; Rafelski et al. 201^. 

The extensive photometric and spectroscopic campaigns 
allowed us to identify galaxy counterpart candidates of two 
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strong Mgii absorbers at z = 0.6915 and 1.4288 towards 
GRB 050820A. The most likely candidates have impact parame- 
ters between 36.8 a nd 42.0 k pc, this being in the expected range 
of Mgii absorbers (IChen et al.ll2010l) . Their properties point to 
young star-forming galaxies, ranging in luminosity from 0.2 
to 0.9 L,. 

Finally, we studied the presence of correlated structures as 
a particular class of intervening absorption-line systems, a phe- 
nomenon observed in the field of GRBs 050730 and 070721B. 
We find evidence for three correlated structures in the field of 
GRB 050730 and one in the field of GRB 070721B. These ab- 
sorbers range in redshift from z = 1.774 to 3.096 and are probed 
over transverse distances between 137 and 161 kpc. All of them 
have associated metal absorption-lines. It is unlikely that the ob- 
jects in these correlated structures are gravitationally bound to 
each other, based on simulations (e.g. Ellison et al. 2007). The 
detection is nevertheless intriguing. The average distance be- 
tween quasar pairs is several times larger. Up to now, intervening 
absorption-line systems with coincidences in redshift have been 
found with separations of several hundreds of kpc up to a few 
Mpc (e.g. Francis & Hewett 1993; D'Odorico et al. 2002), but 
only few with a separation of ~ 150 kpc ( Ellison et al.ll2007ii) . 
These serendipitous discoveries are important for studying the 
absorption cross-section of intervening absorption-line systems, 
their correlation length and their implications on galaxy groups 
in the early Universe. 
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6000 6500 
Observed wavelength (A) 

Fig.A.l. Smoothed spectrum of the galaxy B9 in the field 
of view of GRB 050820A. Based on the identification of [On] 
^3727 and [Oiii] ^5007 the redshift is z = 0.428. The galaxy 
is not related to any of the intervening absorbers seen in the af- 
terglow spectrum. The error spectrum is overplotted. Regions of 
strong telluric features (atmosphere transparency < 20%) were 
not used in the spectral analysis and are marked by small boxes 
at the top. 
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Fig.A.2. Galaxy B2 in the field of GRB050820A. Top: 
Spectrum acquired with VLT/F0RS2. The identified Can H&K 
absorption lines and the B aimer break are marked. The inset 
zooms in on the position of the both absorption lines. Based 
on the identification of the three absorption features the red- 
shift of the galaxy is z = 0.693. The error spectrum is over- 
plotted. Regions of strong telluric features (atmosphere trans- 
parency < 20%) are not used in the spectral analysis and are 
marked by small boxes at the top. Bottom: SED from g'-band to 
5.8 yum. The observed data points (corrected Galactic extinction) 
are shown as circles with error bars. The solid line displays the 
best fit model of the SED (x^ = 16.3; number of filters = 9). 
The model predicted magnitudes are superposed (squares). The 
gray area highlights the interval that is covered by the spectrum 
above. 
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Fig. A.3. Similar to Fig. IA.21 but for galaxy B3 from g'-band to 
5.8 Jim. The fit quality is;^^ = 7.3 for 9 filters. 
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Fig. A.4. Similar to Fig. IA.2I but for galaxy B5 from g'- to Ks- 
band. The fit quality is;^^ = 1.5 for 5 filters. 
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Fig. A.5. Similar to Fig. IA.2I but for galaxy B6 from g'-band to 
5.8 jjm. The fit quality is;^^ = 1.6 for 9 filters. 
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